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20.  Abstract  (cont’d). 

.  size  distribution,  when  normalized  to  the  grain  diameter,  was  very  similar  for  all  specimens  tested.  The  results  indi¬ 
cate  that  the  average  crack  size  is  approximately  one-half  the  average  grain  diameter  over  the  stated  grain  size 
range.  A  dislocation  pileup  model  is  found  to  adequately  predict  the  onset  of  internal  cracking.  The  work  em¬ 
ployed  acoustic  emission  techniques  to  monitor  the  fracturing  activity.  This  information  shed  light  on  the  time 
and  strain  at  which  the  fracturing  began  and  when  the  peak  fracturing  rate  occurred.  Other  topics  covered  in  this 
report  include  creep  behavior,  crack  healing,  the  effect  of  stress  level  on  fracture  size  and  the  orientation  of  cracked 
grains.  Theoretical  aspects  of  the  grain  size  effect  on  material  behavior  are  also  given. 
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EFFECT  OF  GRAIN  SIZE  ON  THE  INTERNAL 
FRACTURING  OF  POLYCRYSTALLINE  ICE 

David  M.  Cole 


Ice  exhibits  brittle  behavior  at  high  tempera¬ 
tures  under  a  variety  of  loading  conditions.  A  key 
factor  causing  this  brittleness  is  the  lack  of  a  suffi¬ 
cient  number  of  independent  slip  systems;  five  are 
required  to  satisfy  the  von  Mises  criterion  for  an 
arbitrary  change  in  shape.  Slip  is  likely  on  the 
basal  and  prismatic  planes  (Goodman  1977),  but 
these  do  not  provide  the  needed  five  systems.  So, 
given  the  inability  of  the  lattice  to  accommodate 
plastic  deformation  by  slip  alone,  and  given  a 
loading  condition  sufficiently  rapid  to  prevent  dif- 
fusional  mechanisms  from  operating  effectively, 
ice  will  develop  cracks.  The  size  and  extent  of  the 
tracks  at  a  given  temperature  depend  on  the  ap¬ 
plied  stress  and  structural  characteristics  of  the 
polycrystalline  aggregate,  such  as  the  orientation, 
shape  and  size  of  the  grains.  Microstructural  dif¬ 
ferences  such  as  those  between  freshwater  ice  and 
sea  ice  also  influence  cracking  activity. 

The  internal  shear  stress  generated  during  either 
tensile  or  compressive  loading  nucleates  cracks. 
Under  certain  conditions,  these  cracks  propagate 
only  a  short  distance  before  coming  to  rest  within 
the  material.  Given  sufficiently  high  stress  levels, 
the  cracks  thus  nucleated  can  propagate  through 
the  material  to  cause  brittle  fracture. 

The  nuclcation  of  stable,  non-propagating 
cracks  is  of  interest  for  a  number  of  reasons. 
These  cracks  are  the  flaws  that  can  propagate 
under  subsequent  tensile  loading.  In  compression, 


when  the  cracks  do  not  propagate,  they  are  re¬ 
sponsible  for  the  gradual  weakening  of  the  struc¬ 
ture  as  straining  proceeds.  Gold  (1970)  noted  that 
ice  passes  directly  from  primary  to  tertiary  creep 
as  a  result  of  the  structural  damage  caused  by  in¬ 
ternal  cracking. 

This  report  concentrates  on  the  effect  of  grain 
size  on  the  internal  cracking  of  polycrystalline  ice 
with  equiaxed  grains.  Relatively  little  research  has 
been  done  in  this  area,  although  considerable 
work  exists  on  the  cracking  of  columnar  grained 
ice,  and  such  work  is  examined  in  detail.  Addi¬ 
tionally,  some  relevant  contributions  regarding 
cracking  in  materials  other  than  ice  are  covered. 

Initial  discussions  center  on  the  root  of  the 
grain-size  dependency  in  material  behavior.  The 
effect  of  grain  size  on  the  ductile/brittle  nature  of 
deformation  is  then  addressed.  Subsequent  sec¬ 
tions  give  attention  to  crack  nucleation  mechan¬ 
isms  and  the  cracking  activity  observed  in  ice.  The 
final  sections  describe  acoustic  emission  tech¬ 
niques  for  crack  detection  and  their  application  to 
the  field  of  ice  mechanics. 


BACKGROUND 

This  work  primarily  examines  the  dislocation 
pileup  mechanism  for  crack  nucleation.  While  the 
operation  of  a  mechanism  based  on  stress  concen- 


trations  arising  from  grain  anisotropy  is  recog¬ 
nized,  it  is  felt  that  the  pileup  mechanism  will 
dominate  at  the  strain  rates  and  temperatures  in¬ 
vestigated  in  this  work. 

Present  research  in  perspective 

The  role  of  grain  boundaries  in  crack  formation 
has  long  been  recognized.  However,  very  little  has 
been  done  to  quantify  the  relationship  between 
grain  size  and  crack  size,  primarily  due  to  the  dif¬ 
ficulty  in  making  the  appropriate  observations  in 
most  materials. 

A  primary  objective  of  this  research  was  to  de¬ 
velop  a  crack  size/grain  size  relationship  for  ice. 
Because  of  its  optical  properties  and  its  propensity 
to  develop  cracks  under  conditions  of  practical 
concern,  polycrystalline  ice  was  ideally  suited  to 
such  a  study  of  internal  cracking  activity.  The  op¬ 
tical  techniques  employed  allow  the  estimation  of 
the  crack  size  distribution  as  well  as  the  number  of 
cracks  per  unit  volume  in  the  tested  material. 

Another  objective  was  to  demonstrate  the  effect 
of  grain  size  relative  to  the  onset  of  internal  crack¬ 
ing.  As  noted  above,  earlier  work  clearly  demon¬ 
strated  the  influence  of  stress  or  strain  rate  on  the 
tendency  of  ice  to  develop  cracks,  but  the  influ¬ 
ence  of  grain  size  alone  in  this  regard  has  not  been 
clearly  demonstrated.  This  objective  was  accom¬ 
plished  by  monitoring  the  extent  of  cracking  in 
specimens  of  increasing  grain  size  while  such  vari¬ 
ables  as  stress,  temperature  and  the  amount  of 
strain  were  held  constant.  The  conditions  for  the 
onset  of  cracking  were  analyzed  in  terms  of  estab¬ 
lished  crack  nucleation  theory. 

Additionally,  this  research  addressed  several 
peripheral  topics  germane  to  the  experimental 
methods  employed  and  to  the  mechanical  proper¬ 
ties  of  ice  in  general.  These  topics  included  the  ob¬ 
servation  of  microcracks  at  various  times  after 
formation  to  monitor  shape  change  (the  crack 
healing  process),  an  examination  of  the  effect  of 
grain  size  on  creep  behavior,  acoustic  emissions 
activity,  and  observations  on  the  orientation  of 
grains  containing  cracks  and  the  orientation  of  the 
cracks  themselves  relative  to  the  axis  of  compres¬ 
sive  stress. 

The  information  obtained  by  the  accomplish¬ 
ment  of  these  objectives  will  be  useful  in  several 
respects.  The  crack  size/grain  size  relationships 
will  enhance  our  understanding  of  the  effect  of 
grain  size  on  the  fracture  strength  of  unflawed  ice. 
Knowledge  of  the  size  distribution  and  number  of 
cracks  will  allow  a  more  precise  examination  of 
the  effects  of  stress/strain  history  on  the  mechani¬ 


cal  properties  of  the  material.  Verification  of  the 
crack  nucleation  model  will  allow  its  application 
with  greater  confidence. 

The  crack  healing  observations  will  be  useful  in 
that  the  results  indicate  the  change  in  crack  geom¬ 
etry  with  time.  This  makes  it  possible  to  assess  the 
period  for  which  a  newly  formed  crack  is  a  signifi¬ 
cant  source  of  internal  stress  concentration. 

Explanations  of  the  grain-size  dependency 

Armstrong  (1979)  pointed  out  the  broad  applic¬ 
ability  of  the  Hall-Petch  relationship  between 
strength  and  grain  size: 

a  =  oo+Kd (1) 

where  a  =  stress 

o9  =  frictional  stress 
K  -  Hall-Petch  slope 
d  =  grain  diameter. 

He  summarized  results  demonstrating  the  validity 
of  the  Hall-Petch  relationship  for  tensile  yield  and 
brittle  fracture  stresses  and  for  the  flow  stress  at 
various  levels  of  strain.  Equation  1  is  essentially 
an  empirical  relationship  and  much  work  has  been 
carried  out  in  efforts  to  develop  a  firm  theoretical 
basis  for  its  veracity. 

Li  and  Chou  (1970)  review  the  major  theoretical 
arguments  that  have  been  put  forth  to  explain  the 
d "  dependency.  Early  work  (see  Stroh  1957  for  a 
useful  summary)  led  to  the  wide  acceptance  of  a 
dislocation  pileup  model  to  explain  the  observed 
grain  size  dependency.  Direct  observations  of  dis¬ 
location  pileups  at  grain  boundaries  made  a  very 
convincing  case  for  this.  Arguments  for  the  dislo¬ 
cation  pileup  model  are  based  on  the  supposition 
that  shear  deformation  passes  from  grain  to  grain 
when  dislocations,  acting  under  an  imposed  stress, 
pile  up  at  a  grain  boundary  and  produce  a  stress 
concentration  that  is  capable  of  producing  slip  in 
the  adjacent  grain. 

A  pileup  at  the  edge  of  one  grain  of  diameter  d 
induces  a  shear  stress  t  at  a  distance  r  in  the  adja¬ 
cent  grain  according  to  the  relationship 

r  =  r,(d/r)v‘  (2) 

where  r,  is  the  applied  shear  stress.  Given  that  a 
stress  is  required  to  generate  slip  in  the  adjacent 
grain,  and  that  a  frictional  stress  r,  must  be  over¬ 
come,  eq  2  may  be  rewritten: 

Tj  =  (rf-  r,)(cf/r)w.  (3) 
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Solving  for  the  applied  shear  stress, 

r.  =  T,  +  ri(r/d)v'.  (4) 

Thus  arises  the  d  *  dependency  according  to  the 
dislocation  pileup  mechanism. 

Interestingly,  as  Li  and  Chou  (1970)  pointed 
out,  materials  in  which  no  pileups  are  observed 
have  been  found  to  obey  the  Hall-Petch  relation¬ 
ship.  This  has  led  to  a  search  for  alternative  ex¬ 
planations  of  the  observed  stress/grain-size  rela¬ 
tionship,  namely  work  hardening  and  grain  boun¬ 
dary  source  theories. 

The  work  hardening  theory  derives  arf'“  de¬ 
pendency  by  using  the  experimentally  established 
fact  that  the  yield  or  flow  stress  is  a  function  of  the 
square  root  of  the  dislocation  density,  q: 

a  =  <Jo  +  Ctflb>J^Q  (5) 

whjre  a  =  a  numerical  constant 

o„  =  the  ordinate  intercept  in  a  plot  of  a  vs 
cf* 

n  =  the  shear  modulus 
b  =  the  Burgers  vector. 

Other  experimental  observations  indicate  that  the 
dislocation  density  at  yield  varies  inversely  with 
grain  size,  thus  explaining  the  d*  dependency. 

The  grain  boundary  source  theory  considers 
grain  boundaries  capable  of  generating  disloca¬ 
tions.  The  length  of  the  dislocation  lines  generated 
in  this  manner  is  directly  proportional  to  the  grain 
boundary  area.  When  this  is  normalized  to  grain 
volume  to  give  a  dislocation  density,  a  ct'  depen¬ 
dency  arises.  Substitution  into  eq  5  again  yields 
the  <t *  dependency. 

Stroh’s  (1957)  work  developed  a  crack  nuclea- 
tion  model  based  on  the  dislocation  pileup  mecha¬ 
nism.  In  order  to  proceed  with  complete  confi¬ 
dence  in  the  use  of  such  a  model,  direct  evidence 
of  dislocation  pileups  in  the  material  in  question 
would  be  necessary.  Sinha  (1978)  presented  photo¬ 
graphic  evidence  of  dislocations  in  polycrystalline 
ice.  Using  an  etching  and  replication  technique,  he 
demonstrated  the  existence  of  dislocation  pileups 
at  grain  boundaries  through  the  observation  of 
etch-pits  on  carefully  prepared  surfaces.  Further¬ 
more,  Sinha’s  results  clearly  indicated  the  glide  of 
basal  dislocations  under  an  applied  stress.  He 
noted  the  appearance  of  dislocations  on  the  f  1 120) 
surface  parallel  to  the  basal  plane.  However,  work 
by  Gold  (1972)  indicates  that  the  pileup  mecha¬ 


nism  may  not  be  the  only  one  to  cause  cracking  in 
ice. 

Gold  (1972)  demonstrated  that  two  independent 
crack  distributions  exist  in  columnar-grained  ice. 
One  was  strain-dependent  and  was  consistent  with 
a  dislocation  pileup  mechanism.  The  other  ap¬ 
peared  to  be  essentially  independent  of  strain,  was 
mainly  composed  of  grain  boundary  cracks  and 
represented  approximately  24%  of  the  total  crack 
population.  Gold  (1972)  speculated  that  cracking 
represented  by  the  latter  distribution  dominated  at 
high  rates  of  loading,  thus  associating  grain  boun¬ 
dary  cracking  with  brittle  behavior.  Furthermore, 
he  suggested  that  the  balance  between  these  two 
crack  distributions  determines  the  transition  from 
ductile  to  brittle  behavior  in  compression. 

The  mechanism  of  this  strain-independent  crack 
distribution  is  not  clear  but  appears  to  be  more 
closely  associated  with  elastic  behavior  than  with 
piastic  behavior.  If  this  is  indeed  the  case,  the  pile¬ 
up  mechanism  should  be  adequate  when  signifi¬ 
cant  plastic  flow  occurs.  However,  its  applicability 
is  liable  to  diminish  as  behavior  becomes  more 
brittle. 

It  is  very  difficult  to  discern  the  crack  nuclea- 
tion  mechanism  from  gross  specimen  observa¬ 
tions.  As  Stroh  (1957)  pointed  out,  the  dislocation 
pileup  model  predicts  the  likelihood  of  cracking  at 
strains  on  the  order  of  those  expected  for  cracks 
caused  by  elastically  generated  stress  concentra¬ 
tions. 

In  light  of  the  above,  while  the  dislocation  pile¬ 
up  mechanism  may  not  be  the  only  source  of  stress 
concentrations  of  sufficient  magnitude  to  generate 
cracks,  it  reflects  the  bulk  of  the  cracking  activity 
of  ice  when  the  behavior  is  not  purely  brittle. 

Grain  size  effects  on  the 
ductile  to  brittle  transition 

Through  its  influence  over  the  internal  distribu¬ 
tions  of  stress,  grain  size  exerts  a  significant  influ¬ 
ence  over  many  aspects  of  material  behavior. 
Most  germane  to  the  present  work  is  the  influence 
of  grain  size  on  the  ductile/brittle  character  of  de¬ 
formation. 

Armstrong  (1970)  explained  the  effect  of  grain 
size  on  the  ductile-to-brittle  transition  in  mild 
steel.  Due  to  the  thermal  effects  on  the  stress  re¬ 
quired  to  cause  either  yielding  or  fracture,  the  fail¬ 
ure  stress  generally  increases  as  temperature  de¬ 
creases.  At  a  constant  strain  rate,  the  material 
undergoes  a  transition  from  ductile  to  brittle  be¬ 
havior  at  some  temperature  Tc.  An  increase  in 


grain  size  lowers  the  peak  stress  experienced  under 
constant  strain  rate  and  increases  Tc.  The  drop  in 
peak  stress  follows  the  slope  of  the  Hall-Petch  re¬ 
lationship.  The  rise  in  Te  results  from  the  relation¬ 
ship  between  grain  size  and  the  temperature- 
dependent  frictional  stress  term  of  the  Hall-Petch 
relationship,  a,. 

At  constant  temperature  and  strain  rate  a  criti¬ 
cal  grain  size  may  be  determined  above  which  the 
material  is  brittle  and  below  which  the  material  is 
ductile.  Stroh  (1957)  arrived  at  a  relationship  be¬ 
tween  transition  temperature  and  grain  size  by  us¬ 
ing  a  stochastic  method: 

1  /  Tc  =  -  V,  (k/u)  log(d)  +  c  (6) 

where  T,  =  transition  temperature 
k  =  Boltzmann  constant 
c  =  a  constant  independent  of  tempera¬ 
ture  and  strain  rate 
u  =  activation  energy. 

More  recently,  Schulson  (1979)  derived  a  relation¬ 
ship  for  the  tensile  case  between  the  critical  grain 
size  and  material  characteristics  of  the  form 

* .  a. 

where  K,c  is  the  critical  stress  intensity  factor.  This 
expression  stems  from  the  fact  that,  at  some  par¬ 
ticular  grain  size,  both  slip-propagation  controlled 
yield  (ductile  behavior)  and  crack-nudeation  con¬ 
trolled  fracture  (brittle  behavior)  are  equally  like¬ 
ly.  Figure  1  shows  the  stress  vs  grain  size  curves 
for  the  ductile  and  brittle  cases.  The  intersection 
defines  the  critical  grain  size.  In  the  present  work, 
grain  size  varies  about  the  critical  grain  size  and 
the  resulting  material  behavior  changes  in  charac¬ 
ter  accordingly. 

The  relationship  between  grain  size  and  7V  can 
be  seen  in  Figure  1 .  A  lower  Te  results  from  a  high¬ 
er  value  of  a.  in  the  Hall-Petch  expression  describ¬ 
ing  curve  I.  This  has  the  effect  of  raising  curve  la 
in  Figure  1  to  curve  lb  and  thus  shifting  the  in¬ 
tercept  with  curve  2  to  a  lower  grain  size.  The  ex¬ 
pression  for  curve  2  is  much  less  sensitive  to  tem¬ 
perature  variations.  Consequently,  it  does  not 
shift  appreciably  and  the  effect  of  temperature  on 
the  point  of  intersection  is  not  significantly  dimin¬ 
ished. 

An  increase  in  grain  size  over  the  critical  value 
brings  about  the  reduction  in  overall  specimen 
strain  prior  to  fracture  often  associated  with  in¬ 
creased  brittleness.  Results  given  by  Mendiratta  et 
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Figure  l.  Stress  /grain  size  relationship  showing 
transition  grain  size  for  ductile  (curve  2)  to  brittle 
(curves  la  and  lb)  behavior.  Shift  from  curve  la  to 
lb  shows  effect  of  decreasing  temperature  on  a,  and  on 
the  critical  grain  size. 


al.  (1976)  show  a  reduction  in  strain  to  fracture  in 
a  titanium  alloy  from  0.21  to  0.02  arising  only 
from  an  increase  in  grain  size.  For  this  change  to 
take  place,  grain  size  was  increased  an  order  of 
magnitude  from  9  to  90  #tm,  and  the  fracture  mode 
changed  from  ductile  dimple  to  brittle  cleavage. 

According  to  work  by  Terlinde  and  Luetjering 
(1982)  grain  size  exerted  an  influence  on  fracture 
strain  of  the  form 

efad-'.  (8) 

In  this  work,  as  in  the  abovementioned  results  of 
Mendiratta  et  al.  (1976),  a  reduction  in  grain  size 
changed  the  behavior  of  a  titanium-aluminum  al¬ 
loy  from  primarily  brittle  to  primarily  ductile, 
with  a  significant  increase  in  failure  strain. 

Nudeation  mechanisms  and  modeling 

Most  current  thought  on  crack  nucleation  stems 
from  a  model  given  by  Zener  (1948).  According  to 
this  model  a  crack  nucleates  when  the  normal 
stress  generated  by  a  dislocation  pileup  reaches  a 
critical  level;  this  causes  the  material  to  fracture 
and  allows  the  dislocations  to  coalesce,  relieving 
the  local  strain  energy.  A  relatively  strong  barrier 
must  be  present  in  order  for  the  pileup  to  build  to 
a  sufficient  stress  to  nucleate  the  crack.  Lattice 
orientation  changes  at  grain  boundaries  or  hard 
inclusions  may  serve  as  effective  barriers. 

Stroh  (1957)  presented  an  extensive  analysis  of 
the  stresses  required  to  nucleate  a  crack.  Stroh 
based  his  development  on  the  concept  of  a  disloca- 
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tion  pileup  on  a  slip  plane,  acted  upon  by  a  shear 
stress,  which  generates  a  sufficient  normal  stress 
in  a  neighboring  grain  to  produce  a  cleavage  frac¬ 
ture.  He  derived  the  expression 

oi  =  3  ir-yAt/8  (1  -v)(  (9) 

for  the  resolved  shear  stress  on  the  slip  plane  and, 

by  using 

f  =  nbn/ir  (1  —  v)a,  (10) 

showed  the  nucleation  condition  to  be 


nocb  =  V,  T2y 


(11) 


where  oc 

y 

V 

( 

b 

n 


resolved  effective  shear  stress  on  the 

slip  plane 

surface  energy 

shear  modulus 

Poisson’s  ratio 

length  of  pileup 

the  Burgers  vector 

number  of  dislocations  in  the  pileup. 


In  eq  9,  coefficients  have  been  determined  for  the 
case  of  the  crack  forming  at  an  orientation  to  the 
slip  plane  which  maximizes  the  stress  on  the  form¬ 
ing  crack.  Stroh  determined  this  angle  to  be  70.5 
He  also  points  out  that  a  crack  length  term  does 
not  appear  in  this  expression. 

In  a  later  work.  Smith  and  Barnby  (1967)  re¬ 
formulated  Stroh’s  approach  to  account  for  the 
effect  of  shear  stress  on  the  nucleation  process  and 
developed  orientation  factors  to  account  for 
geometries  other  than  Stroh’s  case  of  maximum 
normal  stress. 

Smith  and  Barnby  (1967)  give  the  nucleation 
condition  for  a  pileup  of  edge  dislocations  of  a 
single  sign  as 


where  F[<t>)  =  (5  +  2cos <t>  -  3cos20)/4  and  the 
corresponding  number  of  dislocations  required 
under  a,  is 


n 


T!y  1 

2afb  F\<t>) 


(13) 


alysis,  the  crack  length  does  not  appear  in  the  nu¬ 
cleation  criterion,  only  the  pileup  length. 

The  fact  that  the  nucleation  condition  does  not 
contain  a  crack  length  term  is  a  key  point.  The 
length  of  the  crack  is  determined  by  both  the  nor¬ 
mal  stress  component  and  the  presence  of  obsta¬ 
cles  to  its  growth  such  as  grain  boundaries.  In 
other  words,  once  the  separation  of  atom  planes  is 
initiated,  it  will  continue  as  long  as  sufficient  nor¬ 
mal  stress  exists  to  propagate  it.  This  would  be  the 
case,  for  example,  in  a  tension  test  if  the  pileup 
were  of  sufficient  size  to  generate  a  Griffith  crack. 
The  background  tensile  stress  could  drive  the 
crack  (nucleated  via  shear  stresses)  through  the 
material  to  cause  fracture.  If  the  nucleated  crack  is 
not  favorably  oriented  to  the  backgound  stress  or 
if  the  stress  is  of  insufficient  magnitude,  it  will 
come  to  rest  within  the  material. 

In  compression,  the  nucleated  cracks  generally 
do  not  propagate.  Initially,  the  background  com¬ 
pressive  stress  generates  shear  stresses  along  favor¬ 
ably  oriented  slip  planes,  giving  rise  to  dislocation 
piieups,  as  in  the  tensile  case.  Once  the  crack  is  nu¬ 
cleated,  the  compressive  stress  is  not  capable  of 
propagating  the  crack.  Instead,  the  crack  comes  to 
rest  when  the  strain  energy  associated  with  the 
pileup  is  dissipated  or  when  the  leading  edge  of  the 
crack  reaches  a  barrier  that  it  cannot  overcome, 
such  as  the  change  in  lattice  orientation  occurring 
at  a  grain  boundary. 

Visual  observations  (St.  Lawrence  and  Cole 
1982,  Currier  1983)  reveal  a  strong  tendency  for 
non-propagating  cracks  to  form  roughly  parallel 
to  the  loading  axis  in  uniaxial  compression  tests  on 
randomly  oriented,  equiaxed  polycrystalline  ice. 
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Nucleation  conditions  for  more  elaborate  con¬ 
figurations  of  dislocation  sign  and  slip  plane-crack 
geometry  are  also  given.  Again,  as  in  Stroh’s  an- 


Figure  2.  Most  favorable  crack 
orientation  (after  Stroh  1957). 


This  is  reasonable  considering  Stroh’s  determina¬ 
tion  of  the  most  favorable  angle  between  the  slip 
plane  and  the  nucleated  crack.  F  igure  2  shows  the 
geometry  of  this  situation.  The  slip  plane  is  taken 
at  an  angle  of  45°  to  the  loading  axis. 

Although  the  slip  plane  (i.e.  basal  plane)  could 
be  at  an  angle  other  than  the  45"  shown,  the 
planes  of  maximum  resolved  shear  stress  will  tend 
to  cluster  about  this  value.  Also,  Smith  and  Barn- 
by  (1967)  have  shown  that,  while  Stroh’s  optimum 
value  of  =  70.5°  is  correct,  o  may  easily  range 
from  0  to  90°  when  shear  stresses  are  considered  in 
the  analysis.  Even  when  these  values  are  used  as  a 
maximum  range  of  crack  orientation,  the  cracks 
thus  nucleated  will  tend  to  lie  within  about  45 "of 
the  stress  axis  and  have  no  strong  tendency  to 
form  perpendicular  to  it  under  uniaxial  stress. 

Characteristic  size  of  nucleated  crack 

In  examining  fracture  mechanisms  in  metals, 
Gandhi  and  Ashby  (1979)  designated  cracking 
with  no  pre-existing  flaw  as  “cleavage  2.’’  Here 
fractures  are  nucleated  by  slip  or  twinning.  They 
noted  that  these  cracks  were  proportional  to  the 
grain  diameter  and  attributed  this  to  control  by 
the  grain  size  of  the  wavelength  of  the  internal 
stress. 

Physically,  this  proportionality  comes  about  as 
a  result  of  the  obstacle  nature  of  the  grain  boun¬ 
dary.  When  a  polycrystalline  aggregate  is  sub¬ 
jected  to,  say,  a  uniaxial  stress,  the  material  expe¬ 
riences  a  uniform  stress  field  in  a  macroscopic- 
sense.  Microscopically,  however,  this  is  far  from 
the  case:  the  internal  stress  and  strain  fields  are 
very  inhomogeneous. 

Irregularities  in  the  stress  and  strain  fields  are 
brought  about,  in  a  pure  polycrystalline  aggre¬ 
gate,  by  crystal  anisotropy  and  by  dislocation 
movement.  Furthermore,  the  internal  stress  field 
is  in  a  constant  state  of  flux  as  highly  localized  de¬ 
formation  accompanies  both  the  buildup  and  dis¬ 
sipation  of  stress  concentrations  within  the  mate¬ 
rial.  The  frequency  with  which  these  stress  concen¬ 
trations  occur  throughout  the  material  depends 
primarily  on  the  size  of  the  constituent  grains  for 
the  following  reasons. 

The  most  likely  site  for  such  stress  concentra¬ 
tions  is  a  grain  boundary  since  it  offers  a  signifi¬ 
cant  obstacle  to  the  propagation  of  shear  defor¬ 
mation  from  grain  to  grain.  The  most  likely  slip 
plane  is  the  uninterrupted  basal  plane  extending 
across  an  individual  grain.  The  slip  plane  length 
may  or  may  not  equal,  but  in  general  will  scale  as, 
the  grain  size.  Thus,  if  deformation  occurs  via  the 


propagation  of  slip  bands,  the  stress  concentra¬ 
tions  and  their  associated  local  stress  fields  occur 
in  the  material  at  spacings  proportional  to  the 
grain  size.  The  slip  occurs  under  the  action  of 
shear  stress.  However,  the  stress  concentration  re¬ 
sulting  from  slip  generates  a  complex  field  of  ten¬ 
sile,  compressive  and  shear  stresses. 

Gandhi  and  Ashby  (1979)  give  the  expression 
for  a  critical  stress  o*  above  which  a  nucleated 
fracture  will  propagate  and  below  which  the  crack 
will  come  to  rest  with  length  proportional  to  the 
grain  diameter,  d: 

n*  =  (EG./xd)  1  (14) 

where  E  is  Young’s  modulus  and  G,  is  toughness. 

This  is  a  propagation  criterion,  not  a  nucleation 
criterion,  in  that  it  assumes  the  nucleation  of  a 
Griffith  crack  proportional  to  cl.  However,  the  use 
of  a  crack  size  on  the  order  of  d  should  be  noted. 

According  to  Stroh  (1957),  a  nucleated  crack 
will  attain  a  length,  when  normal  background 
stresses  arc  absent,  determined  by  the  number  of 
dislocations  which  enter  it.  Once  the  crack  is  nu¬ 
cleated.  dislocations  enter  it  more  easily  because 
the  back  stress  of  the  pileup  is  relieved.  The  more 
dislocations  that  enter,  the  wider  and  hence  the 
longer  the  crack  becomes.  In  the  compression 
case,  the  only  driving  force  for  the  crack  is  the 
rapidly  relaxing  force  from  the  dislocation  pileup. 
Thus,  the  length  of  the  crack  primarily  is  a  func¬ 
tion  of  the  number  of  dislocations  causing  it  to 
nucleate. 

Generally,  the  analytical  approach  has  been  to 
assume  that  the  favorably  oriented  slip  planes  are 
activated  most  frequently,  and  these  will  in  turn 
nucleate  cracks  most  easily.  If  these  slip  planes 
have  a  characteristic  length,  say  on  the  order  of 
the  grain  diameter,  under  a  given  nominal  stress 
they  will  all  tend  to  contain  about  the  same  num¬ 
ber  of  dislocations.  The  associated  cracks  will  thus 
tend  to  have  a  charcteristic  length  (see  Gold  1966, 
for  example).  Stroh  (1957)  has  indicated  that 
cracks  will  nucleate  and  propagate  to  a  length  on 
the  order  of  f,  the  pileup  length,  in  the  absence  of 
other  driving  stresses. 

However,  the  above  should  be  stated  more  pre¬ 
cisely  in  terms  of  distributions  of  the  quantities 
under  consideration  rather  than  average  or  char¬ 
acteristic  values.  Briefly,  an  estimate  of  the  distri¬ 
bution  of  nucleated  crack  sizes  can  be  obtained  if 
an  appropriate  distribution,  rather  than  an  aver¬ 
age  value,  is  used  to  represent  the  slip  plane 
length.  However,  for  the  purpose  of  demonstra- 


tion  of  the  relationship  between  the  dislocation 
pileup  size  and  the  nucleated  crack  size,  average 
quantities  are  used. 

As  mentioned  above,  grain  boundaries  can  limit 
crack  length.  In  the  case  of  a  uniaxial  compression 
test,  crack  orientation  can  result  in  little  or  no  ten¬ 
sile  stress  normal  to  the  crack  face,  and  a  forming 
crack  may  not  have  a  sufficient  driving  force  to 
overcome  the  crystal  reorientation  at  a  grain 
boundary.  Additionally,  Cottrell  (1958)  viewed 
the  grain  boundary  as  a  likely  stopping  point  for  a 
nucleated  crack  because  a  change  in  orientation  of 
the  cleavage  plane  effectively  represents  a  region 
of  higher  surface  energy  to  the  propagating  micro¬ 
crack.  Thus,  since  grain  boundaries  are  both  likely 
nucleation  and  termination  sites  for  cracks,  the 
crack  size  is  expected  to  correlate  with  the  grain  di¬ 
mension. 

Actual  crack  length  distribution  data  are  un¬ 
common  in  the  literature.  However,  work  by 
McMahon  and  Cohen  (1965)  shows  crack  size  bar 
graphs  for  F4  ferrite  after  repeated  straining  in 
tension.  They  found  that  microcracks  approxi¬ 
mately  equal  to  or  less  than  the  grain  diameter 
formed  first,  and  cracks  up  to  three  times  the 
grain  diameter  formed  after  several  loading  cycles. 
Under  certain  test  conditions,  twin  formation  was 
prevalent  and  the  authors  attributed  a  reduction  in 
the  number  of  large  cracks  to  the  obstacle  nature 
of  the  twins.  Interestingly,  rough  calculations 
based  on  the  bar  graphs  of  McMahon  and  Cohen 
(1965)  indicate  that  the  average  crack  size  is  slight¬ 
ly  over  one  grain  diameter  both  with  and  without 
twin  formation.  They  also  note  that  small  cracks 
continue  to  form  when  the  large  cracks  begin  to 
appear  as  the  number  of  stress  cycles  increases. 

Gold  (1966)  develops  a  quantitative  approach  to 
the  relationship  between  grain  size  and  nucleated 
crack  size  by  using  theory  developed  by  Stroh 
(1954)  and  Bullough  (1964).  Gold  performed  uni¬ 
axial  compression  tests  on  columnar  grained  ice 
and  made  detailed  observations  on  the  size  and 
number  of  microcracks  formed  during  testing. 
Gold's  analysis  considers  the  energy  of  a  cracked 
dislocation  under  an  applied  stress,  and  uses  an 
energy  balance  method  that  leads  to  the  determin¬ 
ation  of  a  critical  or  Griffith  crack.  The  concept  of 
a  cracked  dislocation  as  explained  by  Bullough 
(1964)  allows  the  development  of  a  fracture 
criterion  given  a  dislocation  pileup  and  an  associ¬ 
ated  in-plane  crack  under  an  applied  stress.  The 
fracture  criterion  is  then  based  on  this  crack 
achieving  a  critical  length  for  propagation.  Gold 
(1967)  considers  two  cases  when  the  crack  forms  in 


a  grain  adjacent  to  the  grain  containing  the  pileup: 
1)  the  crack  forms  at  an  angle  to  the  slip  plane 
containing  the  pileup  and  2)  the  crack  forms  in 
plane  with  the  pileup. 

Based  on  the  equations  of  Stroh  (1954)  the 
energy  associated  with  case  1  is 

...  n'b'f t  .  4  L  a„nba 
W  =  T  71  ;  In  —  -  — =- 

4t(1  -  u)  a  2 
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where  W  =  energy  of  the  crack  per  unit  length 
L  =  effective  radius  of  influence  of  the 
dislocations  (L  >  a) 

<r„  =  tensile  stress  perpendicular  to  the 
plane  in  which  the  crack  forms 
a  =  crack  width. 


For  case  2  (Bullough  1964), 
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Gold  (1966)  derives  critical  values  of  crack  width  a 
for  the  two  cases: 


Case  1 


_  2ny 
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Case  2 
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After  assuming  that  a„  is  equal  to  but  of  oppo¬ 
site  sign  than  the  applied  axial  compressive  stress. 
Gold  arrives  at  values  of  a„t,  =  5.7  xlO'4  and  11.4 
x  10'*  m  for  cases  1  and  2,  which  are  in  reasonable 
agreement  with  the  experimental  observations. 
Gold  (1967)  performed  tests  on  replicate  speci¬ 
mens  and  thus  did  not  address  the  issue  of  grain 
size  effects  on  the  cracking  activity  of  the  ice.  In¬ 
deed,  if  the  assumption  regarding  the  value  of  <r. 
being  equal  and  opposite  to  the  applied  stress  is 
maintained,  the  above  expressions  for  critical 
crack  size  are  independent  of  grain  size.  This  is  a 
reasonable  result  when  crack  propagation  is  con¬ 
sidered.  The  cracked  dislocation  is  viewed  as  a  pre¬ 
existing  flaw  and  examined  in  terms  of  its  poten¬ 
tial  to  propagate  under  a  given  stress  field.  The 
potential  (or  likelihood)  for  propagation  is  a  func¬ 
tion  of  material  constants  and  flow  characteris- 
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tics,  but  not  specifically  of  grain  size.  Grain  size 
exerts  only  an  indirect  influence  in  that  it  has  an 
effect  on  the  production  of  flaws  in  the  material. 
Thus,  if  consideration  begins  with  the  material  in 
a  flawed  state,  grain  size  is  not  a  primary  consid¬ 
eration. 

The  material  in  the  present  work  is  considered 
to  be  unflawed  and  the  grain  size-dependent  nucie- 
ation  equations  given  earlier  apply.  In  a  subse¬ 
quent  section,  a  simple  method  is  used  to  relate  the 
nucleated  crack  size  to  the  grain  size  through 
strain  energy  and  surface  energy  considerations  in 
a  manner  similar  to  the  above. 

Cracking  in  ice 

A  series  of  papers  by  Gold  (1960;  1965a,b; 
1966;  1967;  1970a,b;  1972;  1977)  represent  the 
most  extensive  investigations  into  the  internal 
cracking  of  ice.  The  experimental  work  primarily 
involves  columnar-grained  ice,  but  many  of  the 
observations  made  are  germane  to  the  behavior  of 
equiaxial-grained  ice. 

In  his  early  work.  Gold  (1960)  noted  the  forma¬ 
tion  of  cracks  parallel  to  the  long  dimension  of 
grains  in  rectangular  ice  specimens.  The  test  mate¬ 
rial  was  grown  to  result  in  random  c-axis  orienta¬ 
tion  in  the  plane  perpendicular  to  the  long  axis  of 
the  grains.  Cracks  formed  parallel  to  the  grain 
boundaries  and  the  planes  of  the  cracks  were  with¬ 
in  45  °  of  the  stress  axis.  A  detailed  analysis  of  a 
number  of  cracks  indicated  that  30%  were  on 
grain  boundaries,  59%  were  transcrystalline  and 
the  remaining  11%  were  of  mixed  character. 

Gold  (1960)  also  noted  a  change  in  the  cracking 
activity  for  stresses  greater  than  approximately  1.5 
MPa  at  a  temperature  of  -10°C.  SpAdmens  were 
tested  under  creep  conditions  and  the  material  was 
columnar-grained  freshwater  ice.  Tne  average 
grain  diameter  perpendicular  to  the  long  axis  of 
the  grains  was  approximately  4  mm.  Below  the 
1.5-MPa  stress,  cracks  were  relatively  sparse  and 
uniformly  distributed;  above  this  stress,  cracking 
activity  increased  significantly  and  the  cracks 
tended  to  cluster  along  planes  of  maximum  shear. 
Cracking  activity  also  tended  to  peak  early  in  the 
tests. 

Continuation  of  work  along  the  same  lines 
(Gold  1967)  demonstrated  that  cracks  generally  in¬ 
volved  only  one  or  two  grains  and  that  they  tended 
to  propagate  either  parallel  or  perpendicular  to  the 
basal  planes.  The  number  of  cracks  that  formed  in 
a  particular  test  depended  mainly  on  stress  and 
creep  strain  levels  and  was  substantially  indepen¬ 
dent  of  temperature.  The  test  material  in  this  work 
was  again  columnar-grained  freshwater  ice  having 


basal  planes  parallel  to  the  long  dimensions  of  the 
grains.  Compressive  loads  were  applied  perpendic¬ 
ular  to  the  long  dimensions  of  the  grains. 

Gold  (1967)  also  suggests  that  the  strength  of  ice 
experiencing  purely  brittle  failure  is  determined  by 
the  level  of  elastic  stresses  that  the  material  can 
sustain.  In  this  case,  rapid  loading  rates  disallow 
significant  plastic  flow.  Plastic  flow  can  occur, 
however,  at  slower  rates  of  loading,  and  internal 
stress  concentrations  capable  of  initiating  cracking 
eventually  develop. 

Gold  (1967)  found  considerable  scatter  in  the 
time  to  first  crack  formation  under  a  given  creep 
stress.  In  general,  the  time  for  a  crack  to  form 
showed  an  exponential  decay  with  increasing  creep 
stress.  Some  straining  occurred  after  load  applica¬ 
tion  during  which  no  cracks  formed.  For  high 
stresses,  some  small  cracks  appeared  upon  and  im¬ 
mediately  following  loading. 

At  a  level  of  creep  strain  between  3  x  10'*  and 
3  x  10'\  large  cracks  (i.e.  greater  than  2  mm  wide 
x  2  mm  long)  began  to  form.  The  rate  of  forma¬ 
tion  built  up  to  a  peak  and  then  gradually  declined 
as  straining  proceeded.  Results  indicated  that  the 
nucleation  of  a  crack  depended  mainly  on  the  level 
of  creep  strain  and  “not  on  factors  controlling  the 
rate  at  which  the  deformation  occurs.”  Gold 
(1967)  also  noted  that  cracks  tended  to  form  in 
grains  having  their  basal  planes  either  perpendic¬ 
ular  or  parallel  to  the  axis  of  applied  stress.  It  was 
in  this  work,  as  mentioned  earlier,  that  Gold  dem¬ 
onstrated  the  applicability  of  a  dislocation  pileup 
mechanism  to  polycrystalline  ice. 

In  subsequent  work,  Gold  (1970a,b)  quantified 
the  cracking  activity  he  observed  in  columnar- 
grained  ice.  He  also  noted,  for  creep  stresses  less 
than  about  1  MPa  at  temperatures  between  -4.8 
and  -31  °C,  that  cracking  was  confined  mainly  to 
primary  creep  and  that  a  clear  secondary  creep 
stage  developed.  For  stresses  over  about  1.2  MPa, 
however,  continuous  cracking  activity  resulted 
and  the  material  passed  directly  from  primary  to 
tertiary  creep  within  2.5  x  I0~’  strain.  He  attrib¬ 
uted  the  onset  of  tertiary  creep  in  his  columnar¬ 
grained  material  to  the  breakdown  of  the  structure 
by  internal  cracking.  However,  it  should  not  be  in¬ 
ferred  from  this  that  cracking  is  necessary  for  ter¬ 
tiary  creep  to  occur  in  ice  in  general.  Mellor  and 
Cole  (1982)  present  test  results  that  show  a  smooth 
transition  from  primary  to  tertiary  creep  in  the 
absence  of  internal  cracking  in  tests  on  equiaxed 
polycrystalline  ice. 

Gold  (1970a,b)  monitored  the  crack  density  by 
counting  the  number  of  cracks  intersecting  a  plane 
perpendicular  to  the  stress  axis.  Values  were  re- 
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ported  in  the  number  of  cracks  per  unit  area.  By 
deforming  specimens  under  various  loads  or  strain 
rates  (as  well  as  at  several  temperatures)  to  given 
levels  of  strain,  Gold  was  able  to  determine  the  in¬ 
crease  in  crack  density  as  straining  proceeded  for  a 
wide  range  of  test  conditions.  These  tests  yielded 
the  following  additional  information.  The  nucle¬ 
ated  cracks  did  not  appear  to  propagate  with  addi¬ 
tional  straining.  The  cracking  rate  depended  on 
stress,  strain  and  temperature.  The  maximum 
cracking  rate  tended  to  occur  between  axial  strains 
of  1 .5  x  10"1  to  2.5  x  10"’.  For  stresses  below  about 
1  MPa,  the  cracking  rate  tended  to  zero  as  strain¬ 
ing  proceeded.  At  greater  stresses,  cracking  con¬ 
tinued  at  a  reduced  rate  after  the  cracking  rate 
maximum  was  reached.  Cracks  were  randomly 
distributed  in  the  ice  at  low  strains  under  all  test 
conditions.  But  at  higher  stresses  they  tended  to 
form  in  bands  or  “fault  planes”  after  the  maxi¬ 
mum  cracking  rate  had  occurred. 

Although  Gold  did  not  observe  fully  brittle  be¬ 
havior  in  these  tests,  he  did  find  a  decrease  in  the 
strain  at  which  the  strain  rate  minimum  occurred 
with  increasing  cracking.  As  cracking  became 
more  severe,  the  strain  associated  with  the  transi¬ 
tion  to  tertiary  creep  decreased  from  levels  over 
101  to  less  than  2.5  xl0‘5.  He  noted,  however, 
that  even  when  the  lowest  strains  were  observed, 
the  material  response  was  still  significantly  ductile 
in  character. 

Additional  work  (Gold  1972)  reinforced  his  ear¬ 
lier  observations  on  cracking  activity.  He  also  de¬ 
veloped  the  stress  dependency  of  cracking  and  in¬ 
vestigated  the  statistics  of  the  cracking  activity.  He 
found  that  the  crack  sites  are  not  truly  random 
throughout  the  specimen,  but  rather  that  the  prob¬ 
ability  of  a  crack  nucleating  in  a  region  decreases 
if  that  region  already  contains  a  crack. 

Using  Weibull  statistics,  Gold  (1972)  inferred 
the  existence  of  two  separate  crack  distributions. 
One,  believed  to  represent  cracks  generated  by  the 
pileup  mechanism,  was  strain  dependent.  The 
other  distribution  represented  cracks  formed  by 
processes  essentially  independent  of  specimen 
strain;  these  cracks  formed  mainly  at  grain  boun¬ 
daries.  The  probability  of  their  occurrence  in¬ 
creased  with  increasing  applied  stress.  Gold  specu¬ 
lates  that,  given  sufficient  stress,  the  type  of 
strain-independent  cracking  could  be  extensive 
enough  to  be  the  sole  cause  of  specimen  failure.  In 
this  work,  Gold  also  noted  that  crack  density  de¬ 
creased  with  temperature  for  a  given  strain  at  con¬ 
stant  stress. 


'VV'  ■"* 

-  r  j  ’  ■  »  V  ’j 


In  a  review  paper.  Gold  (1977)  pointed  out  the 
need  for  an  increased  understanding  of  the  factors 
influencing  the  cracking  activity  in  ice  as  it  relates 
to  the  ductile-to-brittle  transition,  emphasizing  ice 
type,  temperature,  loading  conditions,  grain  size 
and  specimen  size.  Some  Soviet  workers  have  con¬ 
ducted  work  along  a  similar  line  to  that  of  Gold. 
Zaretsky  et  al.  (1976)  presented  the  results  of  a 
study  on  microcrack  formation  in  columnar¬ 
grained  ice.  As  in  Gold’s  work,  load  was  applied 
perpendicular  to  the  long  axes  of  the  grains  and  the 
c-axes  were  randomly  oriented  in  a  plane  perpen¬ 
dicular  to  the  long  axes  of  the  grains.  They  relied 
heavily  on  the  acoustic  emissions  (AE)  monitoring 
technique  to  quantify  the  cracking  activity.  This 
technique  was  first  used  on  ice  by  Gold  (1960), 
who  subsequently  abandoned  it  and  relied  on  vis¬ 
ual  methods  to  estimate  the  number  of  internal 
fractures. 

The  AE  technique  employs  piezoelectric  trans¬ 
ducers  to  monitor  stress  waves  generated  by  the  in¬ 
itiation  of  a  microcrack.  The  intensity  of  the  stress 
wave  is  assumed  proportional  to  the  magnitude  of 
the  event  that  generates  it.  Electronic  devices  an¬ 
alyze  the  transducer  output  and  characterize  the 
signals  in  various  ways,  depending  on  the  level  of 
sophistication  of  the  particular  system.  A  subse¬ 
quent  section  examines  the  AE  method  in  greater 
detail. 

Zaretsky  et  al.  (1976)  assumed  a  one-to-one  cor¬ 
respondence  between  acoustic  pulses  and  crack 
formation.  Furthermore,  the  amplitude  of  the  AE 
pulse  was  taken  as  proportional  to  the  area  of  the 
crack  that  generated  it.  These  assumptions  were 
substantiated  in  subsequent  work  (Zaretsky  et  al. 
1979). 

Experimentally,  the  Soviet  workers  found  much 
the  same  ice  behavior  as  did  Gold.  Zaretsky  et  al. 
(1976)  found  a  threshold  stress  for  crack  nuclea- 
tion  (denoted  as  as),  The  cracks  tended  to  form 
along  the  grain  boundaries  of  the  columnar- 
grained  test  material.  The  ice  deformed  primarily 
in  two  dimensions— as  also  noted  by  Gold.  Com¬ 
plete  breakup  of  the  specimens  occurred  at  some 
appropriate  level  of  crack  density.  The  number 
and  rate  of  formation  of  microcracks  increased 
with  increasing  applied  stress. 

Zaretsky  et  al.  (1976)  developed  an  equation  for 
the  short-term  ice  creep  rate  in  terms  of  the  ac¬ 
cumulated  number  of  defects  (microcracks),  stress 
and  two  parameters.  By  relying  heavily  on  there 
being  a  relatively  constant  number  of  cracks  at  the 
point  of  “breakup,”  they  developed  expressions 
for  the  time  to  break  up  under  a  given  stress. 
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By  coupling  ice  straining  solely  with  the  occur¬ 
rence  of  internal  cracking  (as  detected  by  AE),  this 
work  inherently  recognizes  cracking  as  the  only 
deformational  mechanism. 

Zaretsky  et  al.  (1979)  expanded  on  much  of  the 
work  presented  in  Zaretsky  et  al.  (1976).  The 
threshold  stress  am  was  viewed  as  the  stress  above 
which  the  “progressive  accumulation  of  structural 
defects  occurs.”  Since  the  accurate  assessment  of 
the  extent  of  internal  cracking  was  critical  to  the 
evaluation  of  the  analytical  expressions  of  this 
work,  Zaretsky  et  al.  (1979)  presented  the  results 
of  a  detailed  petrographic  analysis  on  tested  speci¬ 
mens.  The  results  showed  structural  changes  (i.e. 
the  breakup  of  large  grains)  as  straining  pro¬ 
ceeded.  This  gave  an  indication  of  the  extent  of  in¬ 
ternal  cracking  since  it  was  crack  formation  that 
broke  up  the  large  original  grains  into  smaller 
grains. 

Zaretsky  et  al.  (1979)  concluded  that  the  thresh¬ 
old  stress  a  ,  in  uniaxial  compression,  is  inde¬ 
pendent  of  temperature.  Furthermore,  from  meas¬ 
urements  of  crack  areas,  it  appeared  that  the  mean 
crack  size  increased  with  the  extent  of  cracking 
(and  thus  with  creep  stress).  Crack  size  was  given 
in  arbitrary  units,  however,  and  thus  a  direct  com¬ 
parison  between  crack  size  and  the  grain  size 
(which  ranged  from  2  to  12  mm)  is  not  possible. 
The  analytical  result  was  an  expression  for  creep 
strain  as  a  function  of  stress,  temperature  and  a 
cracking-related  term  based  on  AE  data. 

Detection  of  internal  fracturing 
by  acoustic  emission  techniques 

Information  provided  by  acoustic  emissions 
(AE)  monitoring  can  contribute  significantly  to 
the  understanding  of  material  behavior.  Micro- 
fracturing  activity  especially  lends  itself  to  inter¬ 
pretation  by  AE  techniques  because  pressure 
waves  generated  during  fracture  formation  are 
easily  detected.  The  main  concern  in  handling  AE 
data  is  that  of  interpretation— determining  the  ap¬ 
propriate  correlation  between  the  characteristics 
of  the  source  event  and  the  recorded  AE  signal.  In 
the  present  case,  the  source  event  is  the  formation 
of  a  microcrack  within  the  ice  and  the  relevant  AE 
characteristic  is  signal  amplitude. 

Evans  (1979)  gave  a  theoretical  treatment  of 
acoustical  pulses  generated  by  microfracture  in 
brittle  solids.  He  showed  that  the  amplitude  of 
such  a  pulse  is  a  function  of  crack  geometry,  ap¬ 
plied  stress,  material  properties,  and  distance 
from  the  source.  The  fracture  event  generates  a 
“ringing”  or  oscillation  in  the  crystal  lattice  that 
decays  in  time  and  has  a  characteristic  frequency 


spectrum.  Frequency  analysis  is  often  used  to  dif¬ 
ferentiate  source  mechanisms  where  more  than 
one  mechanism  is  operating.  In  the  present  case, 
however,  only  the  dislocation  pileup  mechanism  is 
assumed  to  be  operating  and  hence  a  frequency 
analysis  is  not  deemed  critical  to  the  investigation. 
The  work  concentrates  primarily  on  the  analysis 
of  AE  pulse  amplitude. 

When  the  size  of  the  AE  source  (i.e.  microfrac¬ 
ture)  varies,  the  peak  AE  signal  strength  varies  as 
well,  all  other  factors  being  equal.  Such  factors  as 
crack  orientation  and  distance  from  the  sensor  can 
cause  significant  differences  in  the  signal  recorded 
for  otherwise  identical  sources.  Thus  it  should  be 
kept  in  mind  that  the  distributions  generally  given 
for  AE  amplitude  reflect  not  only  variations  in  the 
source  itself  but,  to  some  extent,  second-order 
variables  as  well. 

AE  amplitude  data  are  best  manipulated  in  the 
form  of  distribution  functions.  Pollock  (1981)  de¬ 
scribed  the  most  commonly  used  distributions  and 
discussed  their  pros  and  cons.  Among  those  pre¬ 
sented  were  the  Weibull  and  the  log-normal  distri¬ 
butions  as  well  as  the  extreme  value  distributions. 
That  work  also  traced  the  development  of  several 
models  developed  specifically  for  AE  data  analy¬ 
sis.  He  emphasized  that  the  amplitude  distribution 
can  be  considered  a  property  of  the  source  mecha¬ 
nism.  This  last  point  is  a  common  thread  in  much 
AE  work. 

Ono  et  al.  (1978)  showed  close  correlation  be¬ 
tween  the  particle  size  distribution  of  Mn-S  inclu¬ 
sions  in  steel  and  the  AE  amplitude  distribution 
found  when  the  particles  fractured  during  testing. 
Thus,  the  particle  size  governed  the  crack  size, 
which  in  turn  governed  the  amplitude  distribution. 
Wadley  et  al.  (1981)  and  Cousland  and  Scala 
(1981)  are  other  examples  of  work  directed  at  link¬ 
ing  acoustic  activity  and  specific  microstructural 
characteristics. 

In  certain  cases,  AE  also  proves  useful  in  eli¬ 
minating  certain  deformational  mechanisms  from 
consideration.  For  example,  in  Cousland  and 
Scala  (1981),  inclusion  fracture  was  observed  in 
tension  testing  and  produced  very  high  amplitude 
emissions.  When  the  same  material  was  tested  in 
compression,  inclusion  fracture  did  not  occur  and 
no  high  amplitude  AE  signals  were  observed,  thus 
giving  a  clear  indication  of  the  source  mechanism 
in  the  tensile  case.  With  such  information  on  the 
material,  the  extent  of  inclusion  fracture  could  be 
reliably  determined  under  other  test  conditions 
without  extensive  metallographic  investigation. 
Note,  however,  that  quantitative  applications  of 
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AE  require  a  correlation  between  specific  charac¬ 
teristics  of  both  the  deformational  process  and  the 
recorded  AE  signals. 

Wadley  et  al.  (1981)  indicated  the  capability  of 
AE  analysis  to  discriminate  between  two  possible 
event  sources  on  the  basis  of  differences  in  the  AE 
signatures.  By  examination  of  cleavage  and  inter¬ 
granular  crack  sizes  they  noted  that  intergranular 
cracks  were  significantly  larger  on  average  and 
they  corresponded  well  with  the  observed  number 
of  high  amplitude  emissions.  The  particular  equip¬ 
ment  settings  used,  however,  prevented  proper  ac¬ 
quisition  of  the  low  amplitude  signals  resulting 
from  the  smaller  cleavage  fractures. 

Several  studies  have  explored  the  grain  size  ef¬ 
fect  on  internal  fracturing  using  AE  techniques. 
Khan  et  al.  (1982)  found  AE  activity  to  increase 
with  grain  size  for  several  types  of  steel.  Scruby  et 
al.  (1981)  found  similar  trends  for  aluminum  and 
an  aluminum-magnesium  alloy. 

The  optical  clarity  of  ice  and  its  propensity  for 
microfracture  under  conditions  of  practical  inter¬ 
est  make  it  ideally  suited  to  study  with  AE  tech¬ 
niques.  Gold  (1960)  recognized  this  fact  and  was 
the  first  to  use  a  piezoelectric  transducer  to  moni¬ 
tor  cracks  in  ice.  Interpretational  difficulties, 
however,  led  him  to  estimate  internal  cracking  by 
direct  visual  means  in  subsequent  work.  The  po¬ 
tential  benefit  of  the  AE  method  was  clear,  but  the 
equipment  of  the  day  did  not  prove  adequate. 

Work  by  Zaretsky  et  al.  (1979),  mentioned  ear¬ 
lier  in  another  connection,  used  AE  data  directly 
in  a  constitutive  relationship  for  ice.  This  was  pos¬ 
sible  because,  for  certain  test  conditions,  the  ac¬ 
cumulated  acoustic  pulses  followed  the  form  of 
the  accumulated  creep  strain.  This  1979  paper  re¬ 
fers  to  Zaretsky  et  al.  (1976)  for  the  development 
of  a  functional  relationship  between  the  AE  sig¬ 
nals  and  the  corresponding  formation  of  micro- 
cracks.  The  expression  for  creep  strain  was  formu¬ 
lated  as  the  product  of  the  number  of  acoustic 
pulses,  the  mean  crack  size  and  a  proportionality 
factor.  The  mean  crack  size  was  determined 
through  an  analysis  of  the  acoustic  event  ampli¬ 
tudes  by  assuming  that  AE  amplitude  is  a  function 
of  microcrack  size. 

Zaretsky  et  al.  (1979)  also  show  a  close  correla¬ 
tion  between  the  microcrack  surface  area  and  the 
number  of  accumulated  defects.  The  expected  sur¬ 
face  area  increased  linearly  with  the  number  of 
acoustic  pulses  recorded.  The  final  expression 
given  for  ice  creep  in  this  work  showed  time- 
dependent  strain  as  a  function  of  defect  accumula¬ 
tion  measured  by  AE,  temperature,  stress  and 
time. 


II 


The  success  of  this  approach  relies  heavily  on 
the  ability  to  determine  precisely  the  number  of 
cracks  occurring  in  time  from  the  AE  data.  This  is 
a  difficult  task  given  the  variability  of  AE  moni¬ 
toring  systems.  Additionally,  this  approach  is 
valid  only  when  processes  other  than  crack  forma¬ 
tion  do  not  significantly  contribute  to  straining. 

More  recently,  St.  Lawrence  and  Cole  (1982) 
and  Cole  and  St.  Lawrence  (1984)  applied  AE 
techniques  to  monitor  microfracturing  in  poly¬ 
crystalline  ice  having  equiaxed  grains  (in  contrast 
to  the  columnar-grained  material  tested  in  the 
abovementioned  works),  initial  grain  size  was  held 
constant  in  these  experiments  at  1.2  mm  as  deter¬ 
mined  by  the  intercept  method.  Equipment  limita¬ 
tions  prevented  a  direct  correlation  between  AE 
amplitude  and  crack  size.  Instead,  AE  activity  re¬ 
corded  at  two  sensitivity  levels  was  only  assumed 
to  parallel  the  actual  cracking  activity.  The  expres¬ 
sion  developed  for  acoustic  activity  showed  a  de¬ 
pendency  on  stress  and  time.  In  the  creep  tests  re¬ 
ported  in  the  former  paper,  the  AE  rate  reached  a 
maximum  at  1.8x10'’  axial  strain  and  then 
dropped  sharply  as  deformation  proceeded.  For 
stresses  of  less  than  about  2.35  MPa,  the  AE  rate 
after  the  initial  4  x  10  ’  strain  was  extremely  low. 
Stresses  over  3.26  MPa,  on  the  other  hand,  in¬ 
duced  considerable  AE  activity  after  the  rate  peak, 
indicating  a  significant  amount  of  additional  mi¬ 
crofracture. 

In  this  study,  the  acoustic  activity  ranged  over 
some  three  orders  of  magnitude  as  stress  increased 
from  0.8  to  3.67  MPa.  Interestingly,  although  the 
test  material  reached  virtually  complete  saturation 
with  internal  cracks,  the  overall  behavior  could  be 
reasonably  described  as  ductile  since  typical  creep 
behavior  was  still  evidenced  and  the  strain  at  the 
creep  rate  minimum  did  not  decrease. 

In  the  constant  rate  of  deformation  tests  report¬ 
ed  in  Cole  and  St.  Lawrence  (1984),  the  highest 
strain  rates  did  bring  about  substantially  brittle 
behavior.  For  stresses  in  excess  of  5  MPa  and 
strain  rates  over  10“  s  ',  characteristic  failure 
strains  dropped  to  values  as  low  as  2.3  x  10  ’. 
Ductile-type  failures  occur  near  10'!  axial  strain. 
Strain  rates  near  10“  s'1  at  -5°C  resulted  in  virtu¬ 
ally  no  visible  cracking.  An  increase  to  roughly 
10“  s  '  results  in  a  significant  loss  in  ductility  as 
indicated  by  the  occurrence  of  both  a  high  degree 
of  internal  fracture  and  a  reduction  in  the  axial 
strain  associated  with  the  peak  stress. 

In  both  the  creep  and  strength  tests  reported, 
the  onset  of  cracking  as  indicated  by  the  AE  activi¬ 
ty  occurred  at  approximately  10"’  axial  strain. 
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Stress  at  the  onset  of  visible  cracking  was  generally 
near  2.0  MPa  at  -5°C. 

In  other  recent  work,  Sinha  (1982)  monitored 
the  acoustic  activity  in  columnar-grained  ice  in 
uniaxial  compressive  strength  tests.  He  noted  fair¬ 
ly  uniformly  distributed  cracks  that  were  compar¬ 
able  in  size  to  the  grains.  Visible  cracking  gener¬ 
ally  began  at  2.4  x  10'4  strain  and  at  stress  near  0.8 
MPa. 

Sinha  (1982)  associated  visible  cracking  with 
acoustic  event  amplitudes  of  79  dB  with  his  partic¬ 
ular  system.  He  used  visual  observation  during 
testing  to  help  establish  this  cutoff  level. 

As  in  St .  Lawrence  and  Cole  ( 1 982),  Sinha  ( 1 982) 
found  some  low  level  AE  activity  at  the  small 
strains  prior  to  the  onset  of  visible  cracking. 


TEST  METHODS 

This  section  describes  the  testing  methods  and 
procedures  employed  in  the  laboratory  work.  The 
specimen  preparation  procedure  and  the  creep 
testing  equipment  have  been  described  in  detail 
elsewhere  (Cole  1979,  Mellor  and  Cole  1982)  and 
are  covered  only  briefly  here.  However,  the  meth¬ 
od  of  grain  size  analysis  and  the  post-test  analysis 
of  internal  cracking  receive  close  scrutiny. 

Specimen  preparation 

The  specimen  preparation  method  developed  by 
Cole  (1979)  produces  polycrystalline  ice  with  ran¬ 
domly  oriented,  equiaxed  grains  and  densities  of 
0.917  ±0.003  Mg/mJ.  Grain  size  can  vary  up  to 
the  practical  limit  established  by  the  mold  size  and 
is  controlled  by  the  grain  size  of  the  seed  crystals. 
The  specimens  are  50.8  mm  in  diameter  and  127 
mm  long. 

The  method  calls  for  filling  a  cylindrical  alumi¬ 
num  mold  with  the  appropriate  size  seed  grains, 
sealing  the  mold  and  applying  a  vacuum  of  13-26 
Pa  for  2.5  hr.  Distilled,  degassed  water  at  0°C  then 
fills  the  mold  under  the  action  of  the  vacuum. 
Once  this  flooding  is  complete,  the  mold  is  placed 
in  a  freezing  coil  which  carries  fluid  from  a 
temperature  bath  at  -5°C.  The  degassed  water  is 
flushed  up  through  the  mold  as  the  radial  freezing 
progresses  at  an  average  rate  of  2.8  /rni/s.  The 
continuous  flushing  helps  prevent  bubble  nuclea- 
tion  and/or  growth  by  keeping  the  dissolved  gas 
concentration  low  in  the  pore  water. 

The  freezing  process  often  results  in  a  thin  col¬ 
umn  of  fine  bubbles  along  the  axis  of  the  speci¬ 
men.  This  occurs  when  the  freezing  process  pre¬ 
maturely  closes  off  the  path  of  the  flushing  water 


Figure  3.  Typical  untested  specimen. 


near  one  end  of  the  specimen.  The  bubbles  form 
when  the  remaining  gas-laden  pore  water  freezes. 

Figure  3  shows  a  typical  fine-grained  specimen 
produced  by  this  method.  The  end  caps  are  fixed 
in  the  mold  to  assure  proper  alignment.  They  are 
made  from  a  fabric-based  phenolic  material.  The 
ice  bonds  well  to  this  material  once  the  factory  fin¬ 
ish  has  been  roughened  to  expose  the  fabric. 

Specimens  emerged  from  this  procedure  near  a 
temperature  of  -5°C  and  were  placed  in  the  creep 
apparatus  at  -5°C  if  they  were  to  be  tested  imme¬ 
diately.  If  short-term  storage  was  required,  they 
were  wrapped  in  several  layers  of  polyethylene 
film  and  placed  in  ice-filled  bags  and  kept  at 
-12°C.  Such  specimens  equilibrated  at  the  -5°C 
test  temperature  for  at  least  24  hr  prior  to  testing. 

Creep  testing  apparatus 

The  creep  apparatus  and  environmental  control 
cabinet  are  described  in  Mellor  and  Cole  (1982). 
The  end  caps  bolt  into  the  base  and  loading  piston 
of  the  test  fixture.  A  pneumatically  actuated  cylin¬ 
der  applies  the  desired  load  to  the  specimen 
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Figure  4.  Creep  testing  apparatus  showing  dis¬ 
placement  transducer  and  mounting  clamps  on 
specimen. 

through  a  50.8-mm-diameter  steel  piston.  The  pis¬ 
ton  is  mounted  in  a  large  linear  ball  bushing  to  en¬ 
sure  virtually  friction-free  movement. 

The  calibration  procedure  associated  the  output 
of  a  transducer,  which  monitored  the  supply  pres¬ 
sure  to  the  actuator,  to  the  load  exerted  by  the  pis¬ 
ton  on  a  standard  load  cell.  This  method  account¬ 
ed  for  all  frictional  losses  in  the  system. 

The  test  fixture  maintained  the  end  caps  parallel 
during  deformation.  Therefore,  only  one  trans¬ 
ducer  was  required  to  monitor  the  axial  deforma¬ 
tion.  A  direct  current  displacement  transducer 
(DCDT)  with  a  linear  range  of  ±3.175  mm  was 
employed.  Two  circumferential  clamps  held  the 
DCDT  core  and  barrel.  Figure  4  shows  the  com¬ 
plete  creep  testing  apparatus  along  with  the  DCDT 
and  mounting  configuration. 

An  analog  to  digital  data  logger  recorded  the 
DCDT  output  along  with  the  output  of  the  pres¬ 
sure  transducer  and  the  time  of  each  reading.  The 
sampling  time  of  the  data  logger  ranged  from  15  s 


at  the  start  of  a  test  to  as  high  as  300  s  at  higher 
strains  and  slow  strain  rates.  A  separate  system 
continuously  monitored  the  test  temperature, 
which  varied  less  than  ±0.1°C  during  testing. 

Crack  length  and  crack  density  measurements 

After  testing,  specimens  were  moved  to  a  -10°C 
work  room  for  sectioning  and  photographing. 
Specimens  were  generally  cut  on  a  band  saw  to 
generate  horizontal  and  vertical  sections  (see  Fig. 
5).  These  thick  sections  were  approximately  10 
mm  thick,  but  thickness  varied  depending  upon 
crack  density.  High  crack  density  required  thinner 
sections  in  order  to  distinguish  individual  cracks. 
Thicker  sections  could  be  used  when  the  crack 
density  was  low. 

The  horizontal  sections,  taken  perpendicular  to 
the  stress  axis,  were  used  to  estimate  crack  densi¬ 
ties  and  to  measure  crack  lengths.  Since  the  cracks 
tended  to  form  parallel  to  the  stress  axis,  the  hori¬ 
zontally  oriented  sections  showed  the  cracks  in  an 
edge-on  view.  From  this  vantage  point,  the  cracks 
generally  appeared  as  well-defined  lines  and  were 
easily  measured.  The  vertically  oriented  sections, 
while  allowing  measurement  of  crack  dimensions 
parallel  to  the  stress  axis  to  a  certain  extent,  did 
not  provide  an  accurate  means  of  counting  and 
measuring  every  crack  in  the  section.  Inaccuracies 
arose  in  this  case  because  some  cracks  were  seen 
face-on  and  tended  to  obscure  the  view  of  cracks 
which  were  located  behind  them  in  the  section. 


Figure  5.  Schematic 
showing  typical  locations 
of  thick  sections  in  the 
cylindrical  ice  specimens. 
Numbers  1-3  are  horizontal 
sections  used  for  crack  den¬ 
sity  measurements. 


The  question  naturally  arises  as  to  whether  the 
measurement  of  crack  length  in  the  horizontal  sec¬ 
tions  is  an  accurate  representation  of  the  true 
crack  length.  Also,  the  validity  of  the  use  of  one 
length  measurement  to  represent  the  size  of  a 
crack  must  be  established.  These  points  will  be  dis¬ 
cussed  in  detail  and  data  will  be  presented  to  show 
the  extent  of  the  error  introduced  by  these  as¬ 
sumptions. 

For  most  specimens,  photographs  of  back-light¬ 
ed  thick  sections  were  taken.  From  these  it  was 
possible  to  count  and  measure  all  visible  cracks  in 
the  section.  The  sections  were  divided  into  roughly 
200-mmJ  sectors  and  each  sector  was  photo¬ 
graphed  with  a  7  x  magnifying  camera.  It  was 
then  possible  to  form  a  mosaic  of  the  section,  and 
from  this  the  number  and  lengths  of  cracks  were 
taken. 

In  some  cases,  when  the  crack  density  was  ex¬ 
tremely  low,  it  was  possible  to  make  direct  meas¬ 
urement  from  the  viewer  of  the  camera,  preclud¬ 
ing  the  need  of  taking  photographs.  Also  in  these 
cases,  a  larger  volume  of  material  was  sampled  be¬ 
cause  it  was  considerably  less  time  consuming  to 
make  the  measurements.  The  volume  of  the  sec¬ 
tion  was  recorded,  and  once  the  cracks  were 
counted,  the  number  of  cracks  per  unit  volume 
was  calculated. 


Crack  healing  measurements 

k  The  thick-sectioning  technique  described  in  the 

J  previous  section  provided  a  means  to  monitor  the 

!  change  in  crack  length  with  time.  After  photo- 

‘  graphing  immediately  after  testing,  two  typical 

thick  sections  were  tightly  wrapped  and  placed  in  a 
|  -5°C  environment,  and  they  were  photographed 

■  several  times  during  a  period  of  nearly  eight 

i.  weeks.  This  was  sufficient  time  to  allow  complete- 

J  ly  isolated  cracks  to  transform  from  their  initial 

■  “penny”  shape  to  oblate  spheroids.  The  lengths 
of  the  cracks  were  taken  from  each  photograph, 
and  special  attention  was  paid  to  the  first  hours  of 
the  healing  process.  The  results  help  to  assess  the 
possible  change  in  crack  length  resulting  from  the 

•  healing  process  that  occurred  between  the  time  of 
the  crack’s  formation  and  the  time  the  length 
measurement  was  made. 

k 

£  Thin  section  photographs 

■,  Photographs  were  taken  of  thin  sections  and 

•  provided  the  means  of  determining  grain  size  and 
of  discovering  any  anomalies  in  the  test  material 
(see  Fig.  6).  Photographs  were  taken  of  both  verti¬ 
cal  and  horizontal  sections  in  some  cases,  al- 

■  though  generally  only  horizontal  sections  were 


taken.  Figure  6  shows  thin-section  photographs 
for  several  specimens  of  various  grain  sizes. 

The  thick  sections  used  for  the  crack  density  an¬ 
alysis  were  trimmed  to  a  suitable  thickness  for  the 
thin-section  photograph  immediately  after  the 
crack  density  measurements  were  taken.  The 
amount  of  time  between  testing  and  the  final  thin 
section  photograph  was  usually  on  the  order  of  2 
to  4  hours.  Significant  grain  growth  was  not  as¬ 
sumed  to  occur  within  this  time. 

Grain  size  determination 

There  are  several  methods  that  can  be  used  to 
estimate  polycrystalline  grain  size.  A  summary  of 
various  methods  is  given  by  Dieter  (1976)  and  they 
are  briefly  described  below. 

Mean  intercept  length 

Grain  diameter  is  found  by  dividing  the  total 
length  of  a  test  line  by  the  number  of  grains  inter¬ 
sected  when  the  line  is  placed  randomly  on  the  sec¬ 
tion.  This  generally  underestimates  the  true  diam¬ 
eter  of  equiaxed  grains,  but  is  accurate  for  colum¬ 
nar  grains  viewed  perpendicular  to  the  long  axes. 

Grains  per  unit  area 

Assuming  constant  size  spherical  grains,  the 
grain  size  may  be  estimated  by 


where  M  is  the  number  of  grains  per  unit  area. 

A  STM  standard  charts 

Grain  size  at  a  fixed  magnification  is  compared 
with  standard  ASTM  grain  size  charts  and  a  grain 
size  number  is  established.  This  method  will  not  be 
considered  in  the  present  work. 

The  apparent  grain  size  in  the  plane  of  the  sec¬ 
tion  can  also  be  estimated  from  measurements  of 
grains  per  unit  area  M.  In  this  case,  we  find  the 
diameter  which  corresponds  to  the  average  area 
per  grain  1/M,. 


Average  area  = 
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This  results  in  a  somewhat  smaller  estimate  of 
grain  size  than  eq  19.  The  test  results  section  gives 
a  comparison  of  the  grain  sizes  obtained  using 
each  of  the  above  methods.  There  are  significant 
differences  in  grain  size  estimates  depending  on 
the  method  used.  Since  the  work  at  hand  requires 
estimates  of  the  true  grain  size,  and  not  merely 
values  that  scale  as  the  grain  size  (such  as  the  re¬ 
sults  of  the  intercept  method),  the  estimates  result¬ 
ing  from  eq  19,  which  give  the  largest  values,  will 
be  used  in  all  analyses.  The  chosen  method  relies 
on  the  assumption  that  the  grains  are  of  uniform 
size  and  spherical  shape.  Neither  of  these  is  true; 
however,  they  appear  useful  because  the  seed 
grains  are  sieved  to  within  ±8°7o  of  the  average 
seed  size  and  the  seed  grains  are  roughly  equiaxed. 

Caution  must  be  exercised  in  comparing  this 
work  with  other  analyses  in  which  grain  sizes  were 
estimated  with  the  intercept  technique.  As  noted, 
the  grain  sizes  calculated  using  eq  19  are  larger 
than  those  found  with  the  intercept  method  for  the 


present  data.  This  increase  is  significant  and 
should  be  taken  into  account  wherever  grain  size 
measurements  are  of  critical  importance.  In  a  re¬ 
lated  area,  it  should  be  mentioned  that  the  method 
used  to  determine  grain  size  will  influence  the 
slope  of  a  Hall-Petch  type  plot. 

Acquisition  of  acoustic  emission  data 

A  microcomputer-based  AE  system  monitored 
the  acoustic  activity  in  all  tests.  The  system  em¬ 
ployed  two  piezoelectric  transducers  mounted  as 
seen  in  Figure  7.  Elastic  bands  attached  to  the 
mounting  shell  hold  the  transducers  in  place.  Ice 
fillets  formed  from  distilled  water  served  to  in¬ 
crease  the  contact  area  between  the  side  of  the 
specimen  and  the  flat  transducer  face.  A  thin  layer 
of  silicone  grease  between  the  transducer  and  the 
ice  assured  good  acoustic  coupling. 

The  AE  system,  a  PA C  3400  by  Physical  Acous¬ 
tics  Corporation,  recorded  characteristics  of  the 
AE  pulses,  but  not  the  actual  pulse  itself.  Figure  8 


Figure  7.  AE  transducer  mounted  on  specimen.  Trans¬ 
ducers  are  placed  on  flat  contact  points. 
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Figure  8.  Idealized  acoustic  emission  waveforms. 


shows  an  idealized  AE  waveform  and  identifies 
the  major  characteristics  recorded  by  the  system. 
The  gain,  or  amplification  level,  and  the  thresh¬ 
old,  or  cut-off  voltage,  together  determine  the 
overall  sensitivity  of  the  system.  For  these  tests, 
the  gain  was  set  at  60  dB,  which  corresponds  to  an 
amplification  of  1000  times  the  signal  sensed  by 
the  transducer.  The  threshold  setting  varied  some¬ 
what  depending  on  the  AE  activity  level.  St.  Law¬ 
rence  and  Cole  (1982)  point  out  that,  in  ice,  both 
visible  cracking  ana  aetormational  processes 
which  result  in  no  visible  discontinuities  generate 


detectable  acoustic  activity.  Higher  amplitude 
events,  however,  are  expected  from  the  visible 
cracks  as  a  result  of  the  greater  strain  energy  asso¬ 
ciated  with  crack  nucleation.  The  settings  used  in 
this  work  were  such  that  the  AE  system  responded 
to  event  amplitudes  somewhat  below  that  resulting 
from  visible  cracks,  thus  assuring  that  all  the  visi¬ 
ble  cracking  events  were  recorded. 

The  AE  amplifier  band-pass  filters  the  signal  in 
the  range  10  to  200  kHz.  Earlier  work  (St.  Law¬ 
rence  and  Cole  1982)  showed  this  range  to  be  suit¬ 
able  for  monitoring  cracks  in  ice. 


Table  1.  Creep  data. 
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•  Grain  size  achieved  by  grain  growth  process. 

NOTE:  Specimens  that  have  no  values  given  for  and  c  at  toy,  were  not  strained  sufficiently  to  ex¬ 
perience  a  strain  rate  minimum. 


PRESENTATION  OF  RESULTS 


Specimen  characteristics 

Table  1  gives  a  list  of  the  specimens  tested,  the 
initial  applied  stress  level,  the  maximum  axial 
strain  before  removal  of  the  load,  and  the  axial 
strain  at  which  the  minimum  strain  rate  occurred. 

Table  2a  shows  the  specimen  grain  sizes  as  de¬ 
termined  by  the  intercept  method  and  two  meth¬ 
ods  based  on  measurements  of  grains  per  unit 
area.  As  noted  earlier,  the  results  given  in  the  third 
column,  found  using  eq  19,  are  used  in  all  subse¬ 
quent  work  to  characterize  the  material.  These 
values  tend  to  be  significantly  larger  (52.5%  on 
average)  than  those  found  with  the  often  used  in¬ 
tercept  method. 

Table  2b  gives  a  comparison  of  the  three  meth¬ 
ods  of  grain  size  estimation  based  on  thin  sections 
of  untested  material.  The  seed  size  refers  to  the 
sieve  size  range  of  the  ice  crystals  used  to  form  the 
specimen.  Note  that  the  intercept  method  yields 
grain  size  estimates  that  are  smaller  than  the  ori¬ 
ginal  seed  grains.  As  discussed  earlier,  the  method 
used  in  this  work  ( d ,)  gives  estimates  that  are 
slightly  larger  than  the  seed  grains,  but  these  esti¬ 
mates  are  reasonable  because  the  average  seed 
grain  diameter  is  expected  to  increase  as  the  grain 
grows  into  the  adjacent  pore  space  during  freez¬ 
ing. 

Microcrack  measurements 

As  described  above,  post-test  observations  yield 
the  number  and  size  of  cracks  in  a  given  volume  of 
material.  When  the  crack  density  was  very  low,  a 
large  volume  of  material  was  sampled,  and  cracks 
were  measured  and  counted  directly  from  the  thin 
section.  Up  to  three  thick  sections  were  evaluated 
from  each  specimen.  These  data  made  it  possible 
to  estimate  the  crack  density  of  the  entire  speci¬ 
men. 

Since  moderate  to  extensive  cracking  levels  re¬ 
quired  photographs  for  accurate  interpretation, 
two  or  three  sections  of  each  specimen  were  pho¬ 
tographed  as  described  earlier.  Generally,  half  a 
longitudinal  thick  section  was  photographed.  As¬ 
suming  radial  symmetry  in  the  crack  distribution, 
these  results  were  used  to  estimate  the  crack  den¬ 
sity  in  the  central  region  of  the  specimen.  The  esti¬ 
mates  of  crack  density  did  not  include  the  ice  near 
the  ends  of  the  specimen  because  the  triaxial  stress 
state  induced  by  the  end  caps  generally  resulted  in 
a  lower  crack  density  in  these  regions.  Thus,  the 
crack  densities  reported  are  assumed  to  be  repre¬ 
sentative  of  the  material  under  a  uniaxial  stress 
state. 


Table  2.  Grain  size  estimates 
and  seed  grain  sizes. 
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3.6 

3.9 

4.8 

58 

4.7 

4.2 

5.2 

59 

4.2 

4.2 

5.5 

60 

1.7 

2.4 

2.9 

61 

1.9 

2.8 

3.4 

62 

3.6 

4.6 

4.7 

63 

2.0 

2.6 

3.2 

64 

2.0 

2.7 

3.3 

65 

1.7 

2.3 

2.8 

69 

1.2 

1.3 

1.5 

70 

1.1 

1.5 

1.8 

71 

3.6 

4.5 

5.5 

72 

3.4 

4.2 

4.8 

73 

4.4 

4.4 

5.4 

74 

1.9 

2.6 

3.2 

75 

2.0 

2.6 

3.2 

76 

2.1 

2.8 

3.5 

77 

2.1 

2.6 

3.2 

78 

2.4 

2.8 

3.5 

79 

3.7 

4.8 

6.0 

b.  Seed  grain  sizes  and  resulting  grain 
size  measurements  on  untested  ice. 


Seed 

grain  size 

d, 

d, 

d, 

(mm) 

(mm) 

(mm) 

(mm) 

2.80-3.35 

2.5 

3.1 

3.8 

4.0-4.76 

3.6 

4.1 

5.0 

•  dr.  intercept  method. 
d average  area  method. 
d,:  uniform  sphere  assumption  (eq  19). 
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Microfractures 


0  2  4  6  8 

Microfracture  Length  (mm) 

Figure  9.  Typical  crack  length  histogram. 


Table  3.  Crack  location. 


Specimen 

d 

(mm) 

%  GB 
cracks 

V,  XT 
cracks 

49 

3.5 

58 

42 

56 

4.2 

41 

59 

65 

2.8 

SI 

43 

71 

5.5 

54 

46 

72 

4.8 

60 

40 

73 

5.4 

47 

53 

Average 

53  ±7 

47  ±7 

Crack  lengths 

Crack  length  measurements  were  a  direct  result 
of  the  post  test  analysis.  Figure  9  shows  a  typical 
crack  length  histogram.  Appendix  A  contains  all 
the  crack  length  data  presented  in  the  form  of  his¬ 
tograms.  The  average  grain  diameter  is  indicated 
for  each  specimen.  In  some  cases  where  a  more 
detailed  examination  of  the  cracking  was  carried 
out,  transgranular  and  grain  boundary  crack 
histograms  are  shown  separately.  Table  3  gives  the 
mean  values  of  grain  boundary  and  transcrystal- 
line  crack  lengths  obtained  from  these  observa¬ 
tions. 

Figure  10  shows  the  average  crack  length  plot¬ 
ted  against  the  grain  diameter  for  all  cracks,  re¬ 
gardless  of  location.  This  Figure  also  shows  plots 
of  the  least-squares  best-fit  curve  for  all  the  points 
shown  along  with  the  theoretical  curve  to  be  dis¬ 
cussed  later.  The  bars  associated  with  each  point 
indicate  a  bandwidth  of  ±  one  standard  devia¬ 
tion. 

Crack  density 

Figure  11  shows  typical  mosaics  of  the  thick- 
section  photographs.  Each  was  formed  from  five 
enlarged  photographs.  Thin  sections  were  general¬ 
ly  used  to  quantify  the  cracking  activity  for  severe¬ 
ly  cracked  specimens  (Fig.  1  la),  since  the  extensive 
network  of  overlapping  cracks  made  interpreta¬ 
tion  of  the  thick  section  photographs  difficult. 


P'X, 


65 

75  X 

47 

r  z' 

; 

77 

£  , 

43X^ 

631  ' 

2  3  4  5 

Mean  Grain  Diameter  (mm) 


Figure  10.  Mean  crack  length  vs  mean  grain  diameter. 


a.  A  highly  cracked  specimen. 


b.  A  moderately  cracked  specimen. 

Figure  II.  Mosaics  formed  from  enlarged  thick  section  photographs.  Cracks 
show  as  black  lines  of  varying  thickness  under  back  lighting. 


The  crack  density  data,  expressed  as  cracks  per  grain.  Note  that  at  a  constant  initial  stress  of  2.0 

unit  volume,  given  in  Table  4  represent  averages  MPa  crack  densities  range  from  zero  to  nearly  one 

of  several  thick  section  observations  in  most  cases.  crack  per  grain  as  grain  size  increases  from  1.5  to 

Table  4  also  gives  crack  densities  in  terms  of  5.7  mm. 

cracks  per  grain.  These  values  come  about  by  di-  Crack  density  also  increases  with  specimen  strain 

viding  the  observed  number  of  cracks  by  the  esti-  for  all  but  the  smallest  grain  sizes.  Figure  13  shows 

mated  number  of  grains  in  the  sections  under  con-  this  dependency  for  several  grain  sizes, 

sideration.  The  calculation  of  the  number  of  In  some  instances  cracking  in  the  large-grained 

grains  is  based  upon  the  grain  sizes  given  in  col-  specimens  was  so  extensive  that  interpretation  us- 

umn  3  of  Table  2a.  Figure  !2a  gives  plots  of  crack  ing  thick  section  photographs  was  very  difficult, 

density,  in  terms  of  cracks  per  unit  volume,  versus  In  these  cases,  thin  sections  were  prepared  and  the 

grain  size  for  several  ax  tl  strain  levels.  Figure  12b  number  of  cracks  per  grain  was  measured  directly, 

gives  the  same  data  plotted  in  terms  of  cracks  per  The  validity  of  this  procedure  was  checked  by 


Table  4.  Results  of  mtcrofracture  observations. 


Mean  Standard 


d  crack  size  deviation  Mean  _ Crack  density 


Specimen 

(mm) 

(mm) 

(mm) 

d 

(cracks/m') 

(cracks/grain) 

69 

1.5 

no  cracks 

_ 

70 

1.8 

0.83 

0.43 

0.46 

6.8x10’ 

2xl0'> 

43 

1.8 

0.95 

0.40 

0.53 

6.84x10’ 

2.09x10'’ 

55 

(1.8) 

1.39 

0.89 

0.77 

4.83x10’ 

1.47x10“ 

47 

2.6 

1.68 

1.30 

0.65 

1.57x10* 

1.44x10“ 

65* 

2.8 

1.95 

1.32 

0.70 

3.90x10’ 

4.48x10“ 

60 

2.9 

1.28 

0.70 

0.44 

1.03x10’ 

0.13 

77 

3.2 

1.92 

0.99 

0.60 

6.16x10* 

0.10 

63 

3.2 

1.29 

0.90 

0.40 

1.53x10* 

2.58x10“ 

74 

3.2 

1.87 

1.28 

0.58 

8.05x10* 

0.27 

75 

3.2 

1.95 

0.96 

0.61 

4.83x10* 

8.53x10“ 

64 

3.3 

1.14 

0.70 

0.35 

2.38x10* 

4.36x10“ 

44 

3.3 

1.77 

1.14 

0.54 

6.71x10* 

0.125 

61 

3.4 

1.26 

0.86 

0.37 

1.76x10* 

3.59x10“ 

76 

3.5 

1.66 

0.89 

0.47 

4.50x10* 

0.10 

78 

3.5 

2.28 

2.94 

0.65 

5.30x10* 

0.12 

49* 

3.5 

1.18 

0.60 

0.34 

1.95x10’ 

0.434 

56* 

4.2 

1.15 

0.61 

0.27 

1.30x10’ 

0.523 

62 

4.7 

2.09 

1.36 

0.44 

4.26x10’ 

2.53 

57 

4.8 

2.85 

2.66 

0.59 

*.  03x10* 

0.23 

72 

4.8 

3.8 

2.05 

0.79 

1.76x10’ 

1.03 

58 

5.2 

2.51 

1.48 

0.48 

4.20x10* 

0.30 

73 

5.4 

2.45 

1.56 

0.45 

1.83x10’ 

0.80 

59 

5.5 

3.85 

1.45 

0.70 

3.18x10’ 

2.27x10“ 

71 

5.5 

3.38 

2.33 

0.61 

1.28x10’ 

1.13 

79 

6.0 

3.83 

1.78 

0.64 

— 

— 

•  Grain  size  achieved  by  grain  growth  process. 


a.  Number  of  cracks  per  cubic  meter. 

Figure  12.  Crack  density  vs  grain  diameter.  Axial  strain 
levels  are  indicated. 


CO 

sO 


b.  Number  of  cracks  per  grain. 

Figure  12  (cont’d).  Crack  density  vs  grain  diameter.  Axial 
strain  levels  are  indicated. 


comparing  the  results  with  those  of  the  usual 
cracks-per-unit-volume  method  discussed  earlier. 
For  specimen  71,  the  thick  section  analysis  gave  a 
crack  density  of  1.13  cracks  per  grain  and  the  thin 
section  analysis  gave  0.92  crack  per  grain,  indicat¬ 
ing  a  reasonable  agreement  for  the  existing  condi¬ 
tions. 

Additionally,  some  observations  were  made 
from  thick  section  photographs  taken  parallel  to 
the  stress  axis.  These  provided  information  on  the 
shape  of  the  cracks  and  their  orientation  to  the 
stress  axis. 

Creep  behavior 

Figure  14  shows  some  representative  strain-time 
plots  for  tests  at  2.0  MPa.  Figures  15a-c  show 
creep  rate  vs  axial  strain  for  all  tests  conducted  at 
2.0-MPa  axial  stress.  Figures  15d-f  show  similar 
plots  of  the  results  of  tests  subjected  to  the  higher 
stress  levels  (2.4, 2.6  and  2.8  MPa).  The  strain  rate 
minima  show  up  clearly  when  the  data  are  plotted 
in  this  manner.  Note  that  the  scale  of  the  strain 


axes  varies  to  accommodate  the  range  of  strain 
found  in  the  different  tests. 

Most  specimens  tested  to  sufficiently  high 
strains  exhibit  typical  creep  behavior.  The  larger 
grain-sized  material  generally  showed  a  rapidly  de¬ 
creasing  primary  creep  rate,  a  brief  minimum  and 
a  tertiary  phase  in  which  the  creep  rate  tended  to  a 
constant  at  higher  strains. 

The  smaller-grained  material  often  showed  a 
brief  period  of  increasing  creep  rate  at  very  low 
strains.  The  primary  creep  rate  reached  a  maxi¬ 
mum  in  the  range  of  10'5  to  2  x  10"’  strain  and  then 
decreased,  developing  a  relatively  broad  minimum 
near  10'1  axial  strain. 

Several  specimens  were  tested  at  stresses  of  2.4, 
2.6  and  2.8  MPa  in  order  to  examine  the  effect  of 
axial  stress  on  the  cracking  activity  over  a  limited 
range.  Creep  data  for  these  tests  are  plotted  in  Fig¬ 
ures  lSd  and  15c.  No  strong  trends  emerge  from 
the  results  of  these  tests.  Results  given  in  Table  4 
show  that  no  significant  changes  occur  in  the  nor¬ 
malized  crack  size,  indicating  that  the  stress  level 
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a.  Specimens  tested  to  large  strains. 


b.  Some  specimens  tested  to  strains  of  0.25  x  10~2  and 
0.5x10 


c.  Specimens  tested  to  strains  of  approximately  10'2. 


Figure  15.  Creep  curves  for  all  tests. 


d.  Some  specimens  strained  to  approximately  l O'1  under  vari¬ 
ous  stresses. 


e.  Additional  specimens  tested  to  strains  of  approximately 
10  '  under  various  stresses. 


Figure  15  (cont'd). 


/.  Specimens  experiencing  grain  growth  prior  to  testing. 


Figure  15  (cont ’d).  Creep  curves  for  all  tests. 


Figure  16.  Minimum  creep  rate  vs  grain  diam¬ 
eter,  o-2  MPa. 


does  not  affect  the  crack  size  -  grain  size  relation¬ 
ship  over  this  range. 

Although  there  is  scatter  in  the  creep  results,  the 
minimum  strain  rate  increases  with  increasing 
stress.  Inspection  of  the  results  given  in  Table  4  in¬ 
dicates  that  the  number  of  cracks  per  unit  volume 
generally  increases  with  stress  as  well.  The  maxi¬ 
mum  fracture  rate  with  respect  to  time  has  a  mild 
tendency  to  occur  at  lower  strains  as  stress  increas¬ 
es. 

Figure  15f  shows  the  creep  curves  of  three  speci¬ 
mens  for  which  a  relatively  large  grain  size  was 
achieved  through  a  grain  growth  process.  These 
specimens  remained  at  a  temperature  of  -2  °C  for 


a  number  of  weeks  and  consequently  experienced 
considerable  grain  growth.  Upon  testing,  they  ex¬ 
hibited  significantly  different  creep  characteristics 
from  the  other  specimens,  developing  no  decel¬ 
erating  primary  creep  phase.  This  phenomenon  is 
discussed  in  the  section  entitled  The  Effect  of 
Grain  Growth. 

Table  1  gives  the  minimum  creep  rates  along 
with  other  information.  The  creep  rate  minima  are 
plotted  in  Figure  16  as  a  function  of  grain  size. 

Table  1  jpves  the  strain  at  which  the  creep  rate 
minimum  e„.  occurred.  Figure  17  shows  the  strain 
at  lm„  as  a  function  of  grain  size  for  these  tests. 
Note  the  sharp  decrease  in  strain  level  with  in¬ 
creasing  grain  size,  indicating  the  loss  of  ductility 
with  increasing  grain  size. 

Crack  healing 

Observations  were  made  on  two  thick  sections 
to  ascertain  the  extent  to  which  crack  healing 
might  influence  the  crack  length  observations.  Fig¬ 
ures  18  and  19  show  two  sequences  of  time  lapse 
photographs  of  the  crack  healing  process.  The  sec¬ 
tions  were  stored  at  -5  °C  and  photographed  peri¬ 
odically.  Figure  18  shows  an  approximately  3-mm 
crack  face  on.  The  time  after  testing  for  each 
frame  is  given  in  the  figure.  Note  the  prominent 
surface  relief  of  the  crack  face.  Instead  of  a  smooth 
planar  surface,  the  crack  face  appears  to  be  com¬ 
posed  of  many  facets  as  indicated  by  the  network 
of  shadow  lines.  These  features  fade  rapidly  as  the 
crack  heals.  The  crack  surface  becomes  smooth 
and  the  void  is  gradually  transformed  into  a  “bub¬ 
ble”  having  the  shape  of  an  oblate  spheroid.  Fig- 
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Figure  17.  Strain  at  minimum  creep  rate  vs  grain  diameter  for  10  specimens. 


ure  19  shows  several  cracks  edge-on.  The  time  se¬ 
quence  is  the  same  as  for  Figure  18.  The  largest 
crack  did  not  heal  appreciably.  The  reason  is  that 
this  void,  unlike  the  others  studied,  was  exposed 
to  the  atmosphere  and  was  thus  filled  with  air.  As 
discussed  in  greater  detail  below,  this  changes  the 
rate  of  the  healing  process  significantly.  Figure  20 
shows  the  measured  crack  lengths  as  a  function  of 
time. 

Slip  plane  length  distribution 

Initially,  several  methods  of  obtaining  a  distri¬ 
bution  for  t,  the  length  of  slip  plane,  based  on 
ideal  grain  geometries  were  considered.  These 
proved  somewhat  unrealistic  in  light  of  the  differ¬ 
ences  between,  say,  an  idealized  circular  or  hexa¬ 
gonal  cross  section  and  the  variety  of  shapes  actu¬ 
ally  observed  in  the  material  (see  thin-section  pho¬ 
tographs,  for  example).  In  order  to  obtain  a  more 
realistic  sampling  of  slip  line  lengths,  the  follow¬ 
ing  method  was  devised.  A  straight  edge  was  placed 
on  a  thin-section  photograph  and  the  distance  be¬ 
tween  the  grain  boundaries  thus  intersected  was 
recorded  as  the  slip  plane  length.  This  was  done 
several  times  on  each  photograph  and  for  photo¬ 
graphs  from  five  specimens  of  varying  grain  size. 
In  all,  419  measurements  were  taken.  The  values 
from  each  specimen  were  normalized  to  the  speci¬ 
men  grain  diameter,  allowing  all  points  to  be 
merged.  Nearly  10%  of  the  values  were  greater 
than  1.0,  indicating  that  several  of  the  potential 
slip  planes  encountered  were  larger  than  the  aver¬ 
age  grain  size.  This  is  to  be  expected  since  the 


grains  are  not  exactly  uniform  in  size.  The  mean 
and  standard  deviation  of  the  normalized  slip 
plane  lengths  are  0.6  and  0.3.  Figure  21  shows  the 
distribution  of  slip  plane  lengths.  Note  that  the 
maximum  length  is  1.5 d.  The  mean  value  of  Q.6d 
is  used  throughout  this  work  when  an  “average” 
slip  plane  length  is  required. 

Acoustic  emission  observations 

Figure  22  shows  typical  AE  results.  Figure  22a 
gives  accumulated  events  versus  axial  strain.  The 
events  are  normalized  to  unit  volume.  The  curve 
exhibits  the  same  shape  seen  in  other  work  (i.e.  St. 
Lawrence  and  Cole  1982)  but  given  the  filtering 
process  used  in  the  present  work,  the  curve  is  in¬ 
dicative  of  the  actual  number  of  fractures  occur¬ 
ring  per  unit  volume.  Figure  22b  shows  the  deriva¬ 
tive  of  the  curve  in  Figure  22a.  Note  that  most  of 
the  fracturing  activity  occurs  at  very  low  strains. 
Table  5  summarizes  the  AE  data;  it  contains  only 
results  from  specimens  for  which  adequate  AE  data 
were  obtained.  In  certain  cases  technical  problems 
occurred  in  the  recording  of  the  AE  data,  which 
prevented  subsequent  analysis.  In  other  cases  ap¬ 
parently  inadequate  specimen-transducer  contact 
caused  the  results  to  be  of  questionable  validity. 

The  onset  of  fracturing,  as  indicated  by  the  AE 
results,  occurred  at  an  average  axial  strain  of 
4.7  x  10  *  with  a  standard  deviation  of  4.3  x  10  *. 
The  maximum  fracturing  rate  with  respect  to 
strain,  [dAE/de]„„,  occurred  at  1.95x10"’  axial 
strain  with  a  standard  deviation  of  9.6x10'*.  The 
maximum  fracturing  rate  with  respect  to  time, 
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c\  240  hours.  d.  1296  hours. 

Figure  18.  Time-lapse  photographs  of  crack  healing,  face  view. 


\dAE/dt)„a,,  occurred  a!  an  axial  strain  of 
1.75  x  10-’  with  a  standard  deviation  of  8.8  x  10'*. 
These  maximum  rate  calculations  unit  specimen 
59,  which  did  not  experience  sufficiently  high 
strains  to  exhibit  a  plausible  maximum.  There  ap¬ 
pears  to  be  no  systematic  relationship  between  the 
strain  at  the  fracturing  onset  and  grain  size. 

Figure  23  shows  the  effect  of  grain  size  on  the 
time  to  the  maximum  fracturing  rates  with  respect 
to  strain  and  time  for  the  2.0-MPa  tests.  The  time 
to  the  rate  maxima  shows  a  strong  dependency  on 
grain  size,  decreasing  roughly  an  order  of  magni¬ 
tude  as  grain  size  increases  from  1.8  to  5.5  mm. 


Table  5  also  gives  the  AE  amplitude  filter 
threshold  and  the  average  amplitude  of  the  filtered 
events.  The  mean  of  the  filtering  threshold  for  all 
the  tests  is  83.7  dB  with  a  standard  deviation  of 
4.7  dB. 

Grain  orientation 

A  number  of  grain  orientation  measurements 
were  made  to  discern  any  trends  in  the  orientation 
of  grains  having  internal  cracks.  Two  thin  sections 
prepared  from  specimen  79  were  examined  in  de¬ 
tail.  The  thin  sections  contained  numerous  grain 
boundary  and  transgranular  cracks.  The  c-axis 
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a.  Accumulated  acoustic  events  per  cubic  meter 
vj  axial  strain. 


b.  Acoustic  events  per  unit  strain  vs  axial  strain. 


Figure  22.  Typical  acoustic  emission  data. 


Table  5.  Results  of  acoustic  emission  observations. 


AE  rale  maximum 


AE  rate  maximum 


Specimen 

AE  filter 

A  verage 

S/d.  dev. 
(dB) 

First  event 

(with  respect  to  strain ) 

(with  respect  to  time) 

threshold 

(dB) 

amplitude 

(dB) 

€ 

( Xio v 

l 

(S) 

[dAE/dcJm<u 

< 

(xlO1) 

/ 

(S) 

[dAE/dt]max 

(s-') 

e 

(xlO7) 

/ 

(S) 

43 

78.5 

82.2 

2.7 

7.07 

16.0 

3.26x10* 

0.252 

330 

1.58x10' 

0.252 

330 

44 

83.9 

87.9 

3.2 

9.6 

15.0 

1.85x10* 

0.329 

195 

2.30x10’ 

0.225 

105 

47 

— 

— 

— 

0.14 

16.7 

6.57x10* 

0.124 

300 

2.82x10' 

0.088 

210 

49 

76.6 

86.2 

4.1 

0.10 

0.6 

4.03  x  10* 

0.240 

60 

— 

— 

S3 

75.5 

83.7 

3.9 

2.69 

15.3 

5.89x10* 

0.237 

420 

2.26x10' 

0.237 

420 

56 

85.7 

87.8 

1.4 

— 

1.0 

1.16x10'° 

0.312 

150 

5.45x10* 

0.312 

150 

58 

88.0 

90.0 

1.3 

3.6 

9.4 

1.39x10’ 

0.104 

51 
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Figure  23.  Acoustic  emission  rate  data. 

orientations  of  selected  clusters  of  grains  were  de¬ 
termined  using  a  universal  stage  according  to 
methods  given  by  Langway  (1958),  and  using  a 
correction  factor  of  1.04  on  the  equatorial  meas¬ 
urements  rather  than  the  tabulated  values  (see 
Kamb  1962).  Appendix  B  gives  the  results  of  these 
measurements. 

The  results  of  measurements  on  all  grains,  both 
those  containing  cracks  and  the  adjacent  uncracked 
grains,  reflect  the  random  orientation  of  the  test 
material.  No  pattern  of  preferred  orientation  of 
these  grains  emerges.  A  fabric  diagram  for  only 
the  grains  containing  cracks  also  yields  no  discern¬ 
ible  pattern. 

An  examination  of  the  relative  orientation  of 
two  grains  having  a  crack  along  their  common 
grain  boundary  revealed  that  the  angle  between 
the  respective  c-axes  ranges  from  33°  to  85°,  with 
a  mean  of  63.3°  and  a  standard  deviation  of  19°. 

The  above  results  are  not  unexpected,  however, 
since  it  is  the  grains  that  contain  the  active  slip  sys¬ 
tems,  and  not  the  cracked  grains,  where  one  would 


expect  to  find  a  preferred  orientation.  As  noted  in 
an  earlier  section,  the  grains  most  likely  to  slip  and 
thus  cause  crack  nucleation  will  tend  to  have  ori¬ 
entations  that  maximize  the  shear  stress  on  their 
basal  planes.  Unfortunately,  it  was  not  possible  to 
determine  which  grains  were  important  in  forming 
the  observed  cracks.  In  all  probability,  the  grains 
that  generated  the  stress  concentration  were  not 
even  in  the  thin  section. 

Theseresultsalsoindicate  that  the  cracks  formed 
on  planes  other  than  the  basal  in  grains  containing 
transgranular  cracks.  Some  grains  were  observed 
to  have  cracks  running  roughly  parallel  to  the 
c-axis. 


ANALYSIS  AND  DISCUSSION 

As  the  previous  section  shows,  the  results  of  the 
testing  program  provide  information  on  a  consid¬ 
erable  range  of  topics  relative  to  the  response  of 
ice  to  uniaxial  compressive  loading.  These  topics, 
which  are  discussed  in  this  section,  may  be  briefly 
summarized  as  follows. 

1)  Crack  length.  The  average  crack  length  scales 
linearly  with  the  average  grain  size. 

2)  Crack  density.  The  number  of  cracks  per  unit 
volume  increases  with  grain  size  for  the  stated  con¬ 
ditions. 

3)  Onset  of  cracking.  The  test  material  showed  a 
transition  from  ductile  flow  with  no  cracking  to 
flow  with  a  considerable  degree  of  cracking  as  the 
grain  size  increased  from  1.5  to  5.7  mm. 

4)  Creep  behavior.  Grain  size  affected  creep  be¬ 
havior  significantly.  As  grain  size  increased,  pri¬ 
mary  and  tertiary  creep  rates  increased  and  the 
strain  at  the  minimum  creep  rate  decreased  signifi¬ 
cantly. 

5)  Crack  healing.  Observations  indicate  that  iso¬ 
lated,  vapor-filled  cracks  undergo  a  fairly  rapid 
healing  process  which  eventually  transforms  them 
into  a  bubble-like  cavity.  Air-filled  cracks  undergo 
healing  at  a  significantly  slower  rate. 

6)  Stress  effects.  Test  results  at  somewhat  higher 
stresses  indicate  that  the  crack  density  increased 
with  stress,  but  that  the  size  of  the  forming  cracks 
was  not  affected  by  the  stress  level. 

7)  Crack  location  and  grain  orientation.  Cracks 
formed  either  at  grain  boundaries  or  across  grains 
with  nearly  equal  probability,  and  the  grains  that 
developed  cracks  exhibited  no  preferred  orienta¬ 
tion  to  the  stress  axis. 

8)  Acoustic  emission  activity.  All  specimens 
emitted  considerable  acoustic  activity.  Visible 
cracks  generated  the  highest  amplitude  events  and 


the  observations  permitted  some  progress  toward 
establishing  an  event  amplitude  level  associated 
with  the  formation  of  a  visible  microcrack. 

This  section  deals  with  the  above  points  and  ad¬ 
dresses  the  effect  of  the  measurement  techniques 
on  the  results  as  well  as  the  level  of  uncertainty  in 
the  results. 

Thick  section  observations 

The  estimation  of  the  size  and  number  of  cracks 
in  the  thick  sections  proved  to  be  a  very  tedious 
process.  Specimens  with  a  relatively  high  crack 
density  were  difficult  to  analyze  from  the  photo¬ 
graphs  alone  (see,  for  example  Fig.  11a).  In  order 
to  improve  accuracy,  it  was  sometimes  necessary 
to  record  the  number  of  cracks  in  each  sector  of 
the  thick  section  at  the  time  the  photograph  was 
taken,  as  well  as  note  some  general  observations 
about  the  network  of  cracks.  This  additional  in¬ 
formation  greatly  improved  the  reliability  of  the 
measurements. 

Since  an  observer’s  acuity  in  this  type  of  task  is 
likely  to  improve  with  experience,  all  the  crack 


counts  and  measurements  were  carried  out  twice. 
In  some  cases  (most  notably  specimens  49  and  55), 
the  photographs  were  of  relatively  poor  quality 
and  thus  these  results  should  be  treated  with  some 
caution. 

Most  specimens  having  low  or  intermediate 
crack  densities  (less  than  approximately  0.5  crack 
per  grain)  displayed  relatively  unambiguous  net¬ 
works  of  cracks  when  photographed.  The  cracks 
were  isolated  and  were  thus  easily  identified  and 
measured.  The  photograph  in  Figure  lib  shows  a 
typical  sector  from  the  thick  section  of  specimen 
44.  As  noted  earlier,  highly  cracked  specimens 
were  evaluated  by  using  thin  sections.  Additional¬ 
ly,  a  careful  examination  of  the  size  and  number 
of  cracks  which  would  appear  ambiguous  in  the 
thick  section  photographs  provided  information 
which  helped  improve  the  reliability  of  the  results. 

As  noted  above,  thick  sections  taken  parallel  to 
the  axis  of  stress  were  not  suitable  for  crack  densi¬ 
ty  measurements  since  it  was  difficult  to  identify 
individual  cracks  in  that  plane  of  observation.  Fig¬ 
ure  24  shows  a  photograph  of  a  typical  tested  spec- 


Figure  24.  Thick  section  photograph  taken  parallel  to  the  axis  of  applied  stress. 
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and  plane  of  crock 


Figure  25.  Distribution  of  the  angle  be¬ 
tween  the  axis  of  compressive  stress  and 
the  plane  of  the  observed  crack  (40  obser¬ 
vations). 


imen.  It  is  a  thick  section  of  material  near  the  cen¬ 
ter  of  the  specimen  and  it  shows  several  cracks 
face-on  or  nearly  so.  When  viewed  from  this 
angle,  the  cracks  near  the  surface  mask  those 
deeper  in  the  section,  making  an  accurate  count 
nearly  impossible.  However,  this  type  of  photo¬ 
graph  is  useful  in  that  it  shows  the  general  shape 
of  the  cracks  and  it  also  provides  a  means  of  esti¬ 
mating  the  orientation  of  the  plane  in  which  the 
crack  forms  relative  to  the  axis  of  applied  stress. 

Reliable  measurements  of  this  angle  can  be 
made  only  on  cracks  that  appear  edge-on  in  the 
photograph.  Consequently,  relatively  few  meas¬ 
urements  were  available  since  only  several  repre¬ 
sentative  photographs  of  these  vertically  oriented 
thick  sections  were  taken.  Figure  25  shows  the 
results  of  40  observations  in  the  form  of  a  histo¬ 
gram.  Note  that  the  cracks  tend  to  cluster  about 
the  vertical  plane.  The  average  angle  is  23 0  with  a 
standard  deviation  of  17.5°.  Ninety  percent  of  the 
observations  fall  within  45°  of  the  stress  axis. 

Interestingly,  a  few  cracks  formed  at  a  relatively 
large  angle  to  the  stress  axis  and  one  crack  formed 
nearly  perpendicular  to  it. 

The  grain  size  vs  crack  size  relationship 

This  section  examines  the  results  of  the  crack 
length  observations  and  addresses  the  errors  intro¬ 
duced  by  the  observational  method  and  crack 
healing.  Quantification  of  the  sources  of  error 
provides  a  means  of  studying  their  effect  on  the  re¬ 
sults.  Initially,  a  theoretical  grain  size  vs  crack  size 
relationship  is  developed;  later  this  relationship  is 
compared  to  the  experimental  results. 


The  relationship  between  grain  size 
and  nucleated  crack  size 
This  section  develops  a  relationship  between 
grain  size  and  nucleated  crack  size  by  using  the 
concepts  of  elastic  strain  energy  and  surface 
energy.  It  is  similar  to  the  approach  of  Gold  (1966) 
except  that  the  strain  energy  term  considers  the 
pileup  as  a  superdislocation  and  no  flaw  exists  un¬ 
til  the  instant  of  nucleation.  The  method  relies  on 
the  assumption  that  all  the  strain  energy  associ¬ 
ated  with  the  pileup  goes  to  form  new  surface  area 
when  the  crack  forms.  There  is  also  some  inherent 
inaccuracy  because  the  expression  used  for  strain 
energy  was  developed  for  the  isotropic  elastic  case. 
Additionally,  it  is  recognized  that  the  results  of 
such  an  analysis  are  a  strong  function  of  the  values 
assumed  for  various  geometric  parameters  such  as 
obstacle  spacing  and  the  width  of  the  slip  plane. 
The  values  used  here  are  in  agreement  with  those 
found  in  the  literature  (e.g.  Gold  1966)  and  no  at¬ 
tempt  has  been  made  to  obtain  a  range  of  values 
or  a  distribution  to  represent  such  quantities,  with 
the  exception  of  the  slip  plane  length. 

The  calculations  are  based  on  a  grain  diameter 
of  2  mm.  A  slip  plane  length  of  (  =  0.6 d  =  1.2 
mm  is  used  as  a  representative  value.  This  comes 
from  using  the  average  value  of  the  normalized 
slip  plane  length  distribution,  which  was  discussed 
in  a  previous  section.  Only  basal  slip  is  considered 
in  this  treatment.  It  is  also  assumed  that  the  dislo¬ 
cation  pileup  extends  across  a  grain  on  the  slip 
plane  so  the  pileup  length  is  equivalent  to  the  slip 
plane  length. 

Using  eq  12  with 

7  =  0.109  J  m'2  (average  surface  ener¬ 
gy  from  Hobbs  1974) 
g  =  3.5  GPa  (Hobbs  1974) 
t  =  1.2x10-’  m 
0  =  70.50° 
v  =  0.33 

the  resulting  nucleation  stress  is 
o.  =  1.05  MPa. 

When  eq  13  is  used  with  b  =  4.52x10  m,  the 
number  of  dislocations  causing  nucleation  is 

n  =  1200. 


Now  the  elastic  strain  energy  of  n  dislocations  in 
a  pileup  is  given  by  Hirth  and  Lothe  (1968)  as 
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where  (  =  length  of  dislocation  pileup  (slip  plane 
length) 

L  =  length  of  dislocation  line 

R  =  half  obstacle  spacing. 

Assuming  that  L  =  d,  i.e.  the  slip  line  extends 
across  the  grain,  that  f  =  0.6 d  (the  average  of  the 
slip  plane  length  distribution)  and  that  R  =  0.5 d, 
the  energy  of  the  pileup  is  found  to  be 

W  =  2.9  xlO'7  J. 

To  a  first  approximation,  if  the  crack  is  assumed 
to  be  circular,  it  will  have  area  A  =  iry/V4  where 
y'  is  the  calculated  crack  size.  The  nucleated  crack 
has  1A  of  new  surface  area.  The  energy  required 
to  generate  this  area  is  7(2/1).  If  the  elastic  energy 
of  the  pileup  is  equated  to  the  energy  required  for 
generating  surface  area, 

fV  =  7(2/1)  (21b) 

it  is  seen  that 

A  =  (W/2y)  =  1.33xlO-‘  mJ 
or 

yj  =  1.30x10-’  m. 

Thus,  a  crack  size  of  1.30  mm  results  on  average 
for  2.0  mm  grain  size  material  when  it  is  subjected 
to  sufficient  stress  to  nucleate  cracks.  The  ratio  of 
crack  size  to  grain  size  in  this  case  is  0.65.  A  plot 
of  this  relationship  (see  Fig.  10)  indicates  the  lin¬ 
ear  relationship  between  grain  size  and  nucleated 
crack  size. 

Since  in  the  above  analysis  the  crack  size  goes  to 
zero  as  grain  size  goes  to  zero,  and  the  ratio  of  0.65 
may  be  interpreted  as  the  slope  of  the  relationship 
between  the  two  variables,  the  theoretical  relation 
between  crack  size  and  grain  size  may  be  written  as 

y.'  =  0.65  d.  (22) 

where  y.'  is  crack  size  (mm)  and  d  -  grain  di¬ 
ameter. 

Note  that  the  above  treatment  assumes  crack 
formation  a!  all  grain  sizes.  Since  this  is  clearly  not 
the  case,  a  crack  nucleation  criterion  must  be  ap¬ 
plied  to  determine  the  threshold  grain  size  for  nu¬ 
cleation,  and  thus  establish  a  lower  limit  to  eq  22. 

Note  that  the  crack  nucleation  theory  examined 
in  an  earlier  section  (eq  12)  provides  a  means  to  es¬ 
tablish  this  threshold  grain  size  for  crack  nuclca- 
tion. 


Observed  crack  length 
and  sources  of  error 

As  noted  earlier.  Figure  10  shows  the  result  of 
the  crack  size  measurements  for  all  tests  along 
with  the  best  fit  line  using  the  linear  regression 
technique.  The  equation  for  this  line  (r1  =  0.75)  is 

y.  =  -  0.44  +  0.67 d.  (23) 

Equation  22,  the  theoretical  prediction,  is  also 
plotted  in  Figure  10.  Note  that  the  test  results  ex¬ 
hibit  the  linearity  predicted  by  the  model.  Interest¬ 
ingly,  the  theory  predicts  nearly  the  same  slope 
found  in  the  regression  equation,  but  the  curve  is 
shifted  to  somewhat  higher  crack  sizes  throughout 
the  range  of  grain  size. 

The  overprediction  of  the  model  probably  stems 
in  part  from  its  simplicity.  It  is  likely  that  energy  is 
exnended  in  overcoming  the  compressive  back¬ 
ground  stress  in  addition  to  creating  new  surface 
area.  This  would  tend  to  decrease  the  resulting 
crack  size  prediction.  Additionally,  failure  of  the 
crack  formation  event  to  dissipate  all  the  available 
strain  energy  could  contribute  to  the  observed  dis¬ 
parity.  In  fact,  this  circumstance  would  help  ex¬ 
plain  the  difference  in  x-intercept  between  the 
theoretical  and  actual  results.  In  this  case,  equa¬ 
tion  21b  would  take  the  form 

W  =  7(2  A)  +  W„  (24) 

where  W„  is  the  residual  strain  energy  after  forma¬ 
tion  of  the  crack.  The  net  effect  would  be  a  de¬ 
crease  in  crack  size  for  a  given  grain  size  and  a 
shift  in  the  x-intercept  from  the  origin  to  some 
small  grain  size. 

Although  this  potential  source  of  error  in  the 
modeling  will  not  be  addressed  further,  two  sourc¬ 
es  of  error  in  the  crack  length  measuring  process 
require  special  attention,  namely  the  effects  of 
crack  healing  and  the  method  of  crack  length 
measurement  on  the  observed  crack  lengths. 

Crack  healing.  When  a  crack  forms,  it  immedi¬ 
ately  begins  to  heal.  This  is  a  result  of  the  thermo¬ 
dynamic  instability  associated  with  the  crack 
geometry.  The  edges  of  the  cracks  have  an  ex¬ 
tremely  small  radius  of  curvature  and  the  crack 
faces  have  a  relatively  large  radius  of  curvature. 
When  the  crack  remains  isolated  from  the  atmos¬ 
phere  after  its  formation,  it  quickly  fills  with 
water  vapor  in  an  attempt  to  reach  equilibrium. 
However,  the  large  differences  in  surface  curva¬ 
ture  within  the  crack  must  be  eliminated  before 
equilibrium  can  be  achieved. 

Equilibrium  results  when  there  is  no  pressure 
difference  between  the  ice  and  vapor  phases  at  all 
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points  on  the  interface.  This  is  clearly  not  the  cas’ 
in  a  newly  formed  crack.  The  variations  in  the  ra¬ 
dius  of  curvature  result  in  variations  in  the  pres¬ 
sure  difference  between  the  two  phases  along  the 
crack  surface  (see  Colbeck  1980).  The  pressure 
difference  variation  drives  a  flow  process  whereby 
material  is  transported  from  regions  of  low  curva¬ 
ture  to  regions  of  high  curvatun  .  This  process 
gradually  brings  the  void  into  a  nominally  spheri¬ 
cal  shape  with  a  relatively  constant  radius  of  curv¬ 
ature.  The  rate  of  the  process  decreases  as  the  dif¬ 
ference  in  curvature  decreases. 

In  the  present  case  of  an  isolated  crack,  the  void 
contains  only  water  vapor.  The  transport  mecha¬ 
nism  is  viscous  flow  in  the  vapor  phase,  which  is 
fast  relative  to  a  diffusion  process.  The  process 
that  the  ice  undergoes  is  sublimation  in  the  true 
sense  of  the  word,  in  that  it  consists  of  both  evap¬ 
oration  from  and  then  condensation  back  to  the 
solid  phase. 

If  the  crack  and  the  surrounding  ice  are  viewed 
as  a  closed  system  with  no  imposed  temperature 
gradient,  heat  must  flow  from  the  surface  of  the 
crack  receiving  material  to  the  surface  losing  ma¬ 
terial  in  order  for  the  sublimation  process  to  pro¬ 
ceed.  Colbeck  (1986)  points  out  that  this  heat  flow 
is  in  fact  the  rate-limiting  factor  in  the  healing  pro¬ 
cess  of  vapor- filled  cracks. 

Now,  if  the  crack  interior  contains  air,  either 
from  being  opened  to  the  atmosphere  or  from  in¬ 
tersecting  a  gas  bubble  upon  nucleation,  the  heal¬ 
ing  rate  is  slowed  considerably  because  the  trans¬ 
port  mechanism  changes  from  viscous  flow  to  dif- 
fusional  flow,  an  inherently  slower  process.  The 
rate  limiting  factor  is  not  heat  flow,  but  rather  the 
diffusional  process  by  which  water  vapor  travels 
from  the  source  surface  to  the  sink  surface.  Thus, 
the  presence  of  another  gas  in  the  void  results  in  a 
significant  retardation  of  the  sublimation  process. 
Consequently,  the  healing  rale  of  a  gas-filled 
crack  is  much  slower  than  that  of  a  crack  contain¬ 
ing  only  water  vapor  (see  Colbeck  1986  for  an  in- 
depth  treatment  of  this  topic).  The  large  crack  in 
the  center  of  the  photographs  in  Figures  19a-d 
was  filled  with  air  during  the  sectioning  process 
and  displayed  virtually  no  healing  during  the  ob¬ 
servation  period  (see  top  curve  in  Fig.  20b). 

The  question  at  hand  is  whether  this  crack-heal¬ 
ing  process  occurs  fast  enough  to  significantly  in¬ 
fluence  the  crack  length  measurements  taken  up  to 
several  hours  after  their  formation.  To  obtain  the 
answer  to  this  question,  a  set  of  measurements 
were  made  as  discussed  earlier  (refer  to  Fig. 
20a, b). 

Since  most  of  the  cracks  observed  in  these  expe¬ 
riments  are  isolated  within  the  specimens,  it  ap¬ 


pears  that  the  fast  healing  rates  associated  with  the 
vapor  diffusion  mechanism  should  be  considered. 
Considering  the  crack  size  vs  time  data,  the  initial 
healing  rate  is  evidently  not  a  strong  function  of 
initial  crack  size.  Over  the  first  24-hr  period  of  ob¬ 
servation,  the  crack  length  reduction  does  not  vary 
systematically  with  the  original  crack  length,  but 
shows  considerable  scatter  for  cracks  of  similar  in¬ 
itial  length.  This  scatter  is  probably  due  to  varia¬ 
tions  in  the  distance  between  faces  (the  width)  of 
the  individual  cracks.  A  “wide”  crack,  i.e.  one 
with  a  relatively  large  distance  between  opposite 
faces,  requires  transport  of  a  greater  volume  of 
material  to  reduce  its  length  a  given  amount  than  a 
narrow  or  sharp  crack.  If  the  vapor  transport 
mechanism  is  assumed  to  operate  at  nominally  the 
same  rate  for  both  wide  and  narrow  cracks,  the 
wider  cracks  would  require  more  time  to  heal  a 
given  amount. 

Since  crack  length  measurements  were  generally 
taken  less  than  4  hr  after  testing,  estimation  of  the 
maximum  amount  of  healing  likely  to  occur  dur¬ 
ing  this  time  interval  is  of  interest.  Some  healing 
could  have  occurred  during  the  tests,  and  thus  the 
average  amount  of  healing  after  5  hr  should  pro¬ 
vide  a  very  conservative  upper  bound  to  the  crack 
length  reduction. 

The  data  indicate  that,  on  the  average,  the 
measured  crack  length  decreases  by  approximately 
10%  of  the  original  length  in  the  first  5  hr  after 
testing.  However,  as  noted  above,  the  absolute  re¬ 
duction  in  length  due  to  healing  appears  to  be  rela¬ 
tively  independent  of  initial  crack  length,  espe¬ 
cially  at  low  elapsed  times.  This  implies  that  the 
crack  lengths  observed  at  a  fixed  time  after  forma¬ 
tion  require  the  addition  of  a  constant  as  the  ap¬ 
propriate  means  to  correct  for  the  effects  of  the 
healing  process.  The  average  amount  of  healing  in 
the  first  5  hr  after  formation  was  0.08  mm  and  this 
value  is  used  as  the  healing  correction  factor. 
Again,  there  are  large  variations  in  the  amount  of 
healing  of  individual  cracks  and  cracks  which  are 
not  isolated  from  the  atmosphere  will  have  a  negli¬ 
gible  healing  rate  over  this  time  interval. 

Since  the  crack  healing  phenomenon  is  a  peri¬ 
pheral  aspect  of  this  work,  the  complexities  of  the 
healing  process  prevent  an  in-depth  examination 
at  the  present  time.  There  are  many  questions  to 
be  answered  in  this  connection:  the  effect  on  heal¬ 
ing  of  crack  location  (i.e.  grain  boundary  or  intra¬ 
crystalline),  the  point  in  the  healing  process  at 
which  the  crack  becomes  insignificant  as  a  cause 
of  stress  concentration,  and  the  manner  in  which 
crack  size  and  shape  affect  the  healing  process. 

Crack  orientation  and  measured  length.  In  gen¬ 
eral,  the  reported  crack  lengths  are  the  longest  di- 
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Figure  26.  Theoretical  prediction  and  observed  relation¬ 
ship  (after  error  corrections )  for  the  average  crack  size/ 
grain  size  relationship. 


mensions  of  the  cracks  projected  on  a  plane  per¬ 
pendicular  to  the  long  axis  of  the  specimen.  Error 
can  result  if  this  projection  does  not  adequately  re¬ 
flect  the  true  maximum  crack  length.  To  assess 
this  error,  60  cracks  were  evaluated  from  thick 
sections  taken  parallel  to  the  long  axis  of  the  speci¬ 
men.  For  each  crack,  the  ratio  of  the  maximum  di¬ 
mension  to  the  dimension  projected  on  a  horizon¬ 
tal  plane  (which  corresponds  to  the  measurements 
taken  in  the  crack  length  analysis)  was  computed. 
The  minimum  value  of  this  ratio  is  one.  The  ratio 
increases  as  the  true  maximum  becomes  larger 
than  the  maximum  projected  length.  The  average 
of  the  ratio  was  1.12  with  a  standard  deviation  of 
0.19,  indicating  that  the  true  crack  lengths  average 
12%  greater  than  the  values  actually  measured. 

Effect  of  error  on  the  observations 

The  influence  of  the  healing  process  and  the 
method  of  crack  length  measurement  may  be  in¬ 
corporated  into  the  results  by  applying  the  appro¬ 
priate  corrections  to  the  regression  equation  (eq 
23)  which  represents  the  test  data. 

Using  the  conservative  estimate  of  5  hr  for  the 
elapsed  time  between  crack  formation  and  meas¬ 
urement,  the  crack  healing  data  indicate  that,  on 
average,  a  crack  heals  0.08  mm  in  this  period.  The 
regression  equation  may  then  be  corrected  by  add¬ 
ing  the  healing  correction  and  multiplying  the  re¬ 
sult  by  the  factor  1.12.  The  factor  1 .  12  corrects  for 
the  12%  difference  between  the  projected  crack 
length  actually  observed  and  the  true  maximum 
length.  Thus 

y?  =  fy, +0.08)  1.12 


where  y*  is  the  corrected  crack  length  and  yc  is  the 
observed  crack  length.  Substituting  the  expression 
for  yr  from  eq  23, 

y*  =  (  -  0.44  +  0.67d  +  0.08)  1.12 


y*  =  -  0.4  +  0.75 d.  (25) 

Figure  26  shows  the  predicted  relationship  (eq  22) 
plotted  along  with  eq  25,  which  is  the  observed  re¬ 
lationship  after  error  corrections. 

Although  the  corrected  curve  has  a  somewhat 
steeper  slope  than  the  theoretical  prediction,  as 
well  as  a  nonzero  y-intercept,  the  agreement  over 
the  indicated  range  in  grain  size  is  very  good. 

As  a  final  note  in  this  regard,  the  crack  length 
histograms  presented  (i.e.  App.  A)  consist  only  of 
the  uncorrected  observations.  The  effect  of  apply¬ 
ing  the  corrections  to  these  data  would  be  a  shift 
to  larger  crack  lengths  and  a  slight  increase  in  the 
range  or  spread  of  the  data. 

Crack  nucleation  condition 

The  results  show  that  under  the  prevailing  test 
conditions,  cracks  begin  to  nucleate  at  grain  sizes 
between  1.5  and  1.8  mm.  Recall  that  the  1.5-mm 
grain-sized  specimen  exhibited  no  cracking  and 
very  slight  cracking  was  observed  in  the  1.8-mm 
grained  specimen  (see  Table  1).  Thus  the  threshold 
grain  size  lor  crack  nucleation  lies  in  this  range  of 
grain  sizes. 
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Given  this  information,  then,  the  validity  of  the 
crack  nucleation  condition  (eq  12)  may  be  tested 
by  solving  for  the  slip  plane  length  and  substitut¬ 
ing  appropriate  values  for  the  various  elements  of 
the  relationship.  If  eq  12  adequately  represents  the 
mechanism  of  crack  nucleation,  the  resulting  criti¬ 
cal  slip  plane  length  for  crack  nucleation  should 
coincide  reasonably  well  with  the  observed  critical 
grain  size  for  crack  nucleation.  The  development 
of  this  correspondence  between  slip  plane  length 
and  grain  size,  as  noted  earlier,  relies  on  the 
assumption  that  shear  strain  propagates  via  slip 
planes  that  extend  completely  across  the  grains. 

The  following  expression  is  obtained  by  solving 
eq  12  for  the  critical  slip  plane  length,  t: 

f  1 

2(1  -v)al  F(4>)  • 

Values  for  y,  n  and  i>  given  earlier  are  used.  How¬ 
ever,  the  values  for  oL,  the  effective  shear  stress  re¬ 
quired  for  crack  nucleation,  and  the  geomet¬ 
rical  parameter,  require  some  discussion. 

Under  the  prevailing  test  conditions,  the  maxi¬ 
mum  resolved  shear  stress  aMKss,  which  occurs  on 
the  most  favorably  oriented  slip  plane,  is  the  axial 
stress  aA  multiplied  by  the  maximum  Schmid  fac¬ 
tor  m  or 

Omrss  —  Qa  •  m. 

Since  oA  =  2.0  MPa  and  the  highest  value  of  m 
is  0.5, 

&MRSS  —  1.0  MPa. 

Now  the  maximum  effective  shear  stress  aM£J,  is 
the  maximum  resolved  shear  stress  less  the  fric¬ 
tional  stress  component  a0,  or 

Omess  =  Omrss  ~  O0. 

a0  represents  the  lattice  resistance  to  dislocation 
movement  in  terms  of  shear  stress.  The  only 
source  for  a  value  of  a0  for  ice  prepared  in  the 
same  manner  as  in  this  work,  and  at  the  same  test 
temperature  of  -5°C,  is  Lim  (1983).  He  reported  a 
value  of  0.56  MPa  for  the  frictional  component  of 
stress,  in  terms  of  axial  stress,  for  polycrystalline 
ice  tested  in  tension  under  an  average  strain  rate  of 
9.4  xlO'7  s'1.  Since  the  present  work  deals  with 
shear  stress,  the  value  of  0.56  MPa  has  been  multi¬ 
plied  by  0.5  and  rounded  to  one  significant  figure, 
yielding  a  value  of  ct0  =  0.3  MPa  for  the  shear 


stress  required  to  overcome  lattice  resistance  to 
dislocation  movement.  The  direction  of  applied 
stress  is  assumed  to  exert  no  influence  on  the  fric¬ 
tional  stress  value. 

Thus,  the  value  of  the  maximum  effective  shear 
stress  ac  to  be  used  in  eq  24  is 

Oe  ~  Omrss  ~  <*o 

=  1.0  MPa  -  0.3  MPa 

or  a£  =  0.7  MPa. 

The  function 

F{<t>)  =  (5  +  2  cos  <t>  -  3cos2<£)/4 

given  by  Smith  and  Barnby  (1967)  provides  a 
means  of  making  the  nucleation  stress  a,  sensitive 
to  the  slip  plane/crack  plane  geometry.  In  his 
treatment,  Stroh  (1957),  as  noted  above,  calculat¬ 
ed  the  optimum  value  of  <t>  to  be  70.5  0  based  on 
the  normal  stress  distribution  associated  with  the 
pileup.  His  results  show  this  point  to  be  relatively 
well  defined.  That  is,  the  probability  of  develop¬ 
ing  a  crack  at  a  larger  or  smaller  angle  <t>  decreases 
sharply  as  the  angle  deviates  from  70.5°. 

By  considering  shear  stress,  Smith  and  Barnby 
(1967)  showed  that,  while  the  optimal  value  of 
70.5°  was  correct,  the  angle  <f>  can  vary  from  0°  to 
somewhat  over  90°  with  relatively  little  change  in 
the  stress  required  for  crack  nucleation.  This  re¬ 
sult  has  the  net  effect  of  making  the  stress  required 
for  nucleation  of  a  crack  less  sensitive  to  local 
geometric  factors. 

Considering  the  likely  range  in  <j>  to  be  0°  to  90°, 
the  minimum  and  maximum  values  of  F(<f>)  are  1 .0 
and  1.34.  These  values  occur  at  <t>  =  0°and  70.5°. 

It  is  now  possible  to  calculate  a  range  in  slip 
plane  lengths  over  which  crack  nucleation  is  possi¬ 
ble  given  the  prevailing  test  conditions.  Equation 
26  yields  the  following  values: 

F(<t>)  =  1.00,  (  =  1.8  mm 

F(0)  =  1.34,  f  =  1.4  mm. 

This  range  in  slip  plane  length  agrees  well  with  the 
range  in  grain  size  (1.5  to  1.8  mm)  over  which  the 
transition  from  no  cracking  to  cracking  was  ob¬ 
served.  Thus,  given  the  assumptions  mentioned 
above,  the  nucleation  criterion  (eq  12)  appears  to 
model  the  observed  ice  behavior  reasonably  well. 

The  fact  that  ice  has  a  relatively  low  coefficient 
of  self-diffusion  £>,  is  undoubtedly  a  major  factor 
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in  the  agreement  between  theory  and  observation 
in  this  case.  D,  for  icc  at  -10°C  (corresponding  to 
a  homologous  temperature  of  0.96)  is  on  the  order 
of  10'"  m!  s'1  (Hobbs  1974)  while  a  typical  value 
for  a  metal  is  on  the  order  of  10  "  m!  s'1 
(Shewmon  1969).  The  lower  capability  for  the  dif¬ 
fusion  of  vacancies  in  ice  enhances  the  material’s 
ability  to  build  and  sustain  the  stress  concentra¬ 
tions  necessary  for  the  nucleation  of  cracks  at  high 
homologous  temperatures.  In  contrast,  at  high  ho¬ 
mologous  temperatures  in  metals,  the  high  degree 
of  vacancy  diffusion  inhibits  pileup  formation  and 
thus  precludes  the  development  of  stress  concen¬ 
trations  necessary  for  crack  nucleation. 

Crack  density  and  specimen  strain 

It  is  useful  to  examine  the  cracking  activity  as  a 
function  of  strain.  The  results  given  in  Figure  13 
are  in  qualitative  agreement  with  the  work  of 
others  who  have  made  direct  observations  of  inter¬ 
nal  cracks  in  ice  (Gold  1970a,  Zaretsky  et  al.  1979). 

Typically,  cracks  begin  to  form  after  a  small 
amount  of  strain  has  occurred.  The  acoustic  emis¬ 
sion  results  given  later  indicate  the  strain  level  for 
the  onset  of  visible  cracking  is  4.7  x  10'*  on  aver¬ 
age.  Interestingly,  this  strain  level  agrees  well  with 
Gold’s  (1970a)  observations  on  columnar-grained 
ice  at  -9.5  °C.  He  found  this  strain  level  to  be  rela¬ 
tively  independent  of  stress  provided  it  was  suffi¬ 
cient  to  cause  cracking. 

The  results  show  that  the  manner  in  which  cracks 
accumulate  with  strain  depends  upon  the  grain 
size  (refer  again  to  Fig.  13).  For  the  small-grained 
material,  most  all  the  cracking  occurs  in  the  first 
10"  strain.  Additional  straining  results  in  a  negli¬ 
gible  increase  in  the  number  of  cracks.  However, 
at  the  intermediate  grain  sizes  (approximately  3 
mm),  the  cracking  continues,  but  at  a  much  reduced 
rate  as  straining  proceeds  beyond  10".  Finally,  the 
largest-grained  material  exhibits  a  very  high  initial 
rate  of  cracking  at  strains  below  10"  and  then  con¬ 
tinues  to  generate  cracks  at  a  significant  rate.  At 
total  strains  of  less  than  3  x  10",  these  specimens 
are  so  completely  saturated  with  cracks  that  they 
appear  opaque.  The  crack  density  can  be  esti¬ 
mated  only  from  thin  sections  for  these  specimens. 

These  results,  brought  about  only  by  an  increase 
in  grain  size,  are  in  good  agreement  with  the 
trends  observed  by  St.  Lawrence  and  Cole  (1982). 
These  resulted  from  an  increase  in  applied  stress, 
while  grain  size  remained  constant  at  1 .2  mm  (esti¬ 
mated  by  the  intercept  method)  and  the  test  tem¬ 
perature  was  -5°C.  In  that  work,  an  increase  in 
stress  brought  about  an  increase  in  cracking  activ¬ 
ity  as  indicated  by  acoustic  emissions  monitoring. 


A  brief  comparison  of  the  influence  of  stress  and 
of  grain  size  on  cracking  shows  that  the  increase  in 
grain  size  in  the  present  work  has  an  effect  on  the 
cracking  activity  that  is  very  similar  to  the  effect 
of  an  increase  in  axial  stress.  This  is  not  surpris¬ 
ing,  however,  since  it  is  well  recognized  that  an  in¬ 
crease  in  either  stress  or  grain  size  generally  re¬ 
duces  ductility. 

Creep  behavior 

Grain  size  effects 

The  character  of  the  creep  behavior  seen  in  Fig¬ 
ures  15a-c  is  similar  to  that  found  in  other  work 
for  the  finer-grained  material  at  similar  stress 
levels  (Mellor  and  Cole  1982,  Jacka  1984). 

The  effect  of  grain  size  on  the  creep  curve  can 
be  seen  by  comparing  specimens  62  and  69  in  Fig¬ 
ure  15a.  An  increase  in  grain  size  from  1.5  to  4.7 
mm  causes  a  drop  in  the  strain  at  the  minimum 
creep  rate  from  10"  to  4.2  x  10".  There  is  a  tend¬ 
ency  for  the  larger-grained  material  to  exhibit 
both  a  faster  drop  in  the  primary  creep  rates  and  a 
faster  increase  in  the  tertiary  rates  than  the  fine¬ 
grained  material. 

Some  of  the  tests  (see  Fig.  15a,  specimens  47,  60 
and  62)  exhibited  a  trend  toward  higher  primary 
creep  rates  with  larger  grain  size.  Duval  and 
LeGac  (1980)  observed  this  trend  in  creep  tests  on 
polycrystalline  ice  at  a  temperature  of  -7°C  and  a 
creep  stress  of  approximately  0.5  MPa.  The  same 
workers  also  noted  that  grain  size  exerted  no  ap¬ 
parent  influence  on  the  “steady  state”  creep  rate 
of  their  test  material.  The  results  of  Duval  and 
LeGac  (1980)  are  at  variance  with  results  reported 
by  Baker  (1978),  who  observed  a  significant  effect 
of  grain  size  on  the  steady-state  creep  rate  of  lab¬ 
oratory-prepared  ice.  Baker  found  the  creep  rates 
exhibited  a  minimum  at  a  grain  size  of  1.0  mm  for 
tests  performed  at  temperatures  of  -7  to  -10°C 
and  at  a  creep  stress  of  approximately  0.56  MPa. 
He  attributed  the  observed  reversal  of  the  influ¬ 
ence  of  grain  size  (from  strengthening  to  weaken¬ 
ing  the  material)  to  a  change  in  the  main  deforma¬ 
tion  mechanisms.  He  reasoned  that  the  grain 
boundary  weakening  resulted  from  diffusional 
processes  operative  for  the  small-grained  material 
and  that  the  strengthening  resulted  from  the  oper¬ 
ation  of  the  dislocation-controlled  creep  mecha¬ 
nism. 

Unpublished  work  referenced  by  Jacka  and 
Maccagnan  (1984)  indicates  no  significant  grain 
size  effect  on  the  minimum  creep  rate  of  lab¬ 
oratory-prepared  ice  over  the  grain  size  range  of 
0.8  to  3.4  mm,  at  temperatures  of  -7  and  -10°C 


I 


B 


40 


and  under  creep  stresses  of  0.3  and  0.26  MPa  octa¬ 
hedral  (0.64  and  0.55  MPa  normal).  These  experi¬ 
mental  conditions  cover  the  range  of  Baker’s 
(1978)  tests,  but  the  grain  size  effect  observed  by 
Baker  is  absent.  Thus,  the  work  by  Baker  (1978) 
remains  unsubstantiated.  The  reason  for  the  dis¬ 
agreement  on  grain  size  effects,  however,  is  not 
clear.  Historically,  a  main  difficulty  in  the  field  of 
ice  mechanics  lies  in  the  effect  of  the  specimen 
preparation  procedure  on  the  mechanical  proper¬ 
ties  of  the  material.  Since  there  is  no  standard 
method  for  specimen  preparation,  workers  gener¬ 
ally  begin  with  a  commonly  used  approach  such  as 
packing  a  mold  with  sieved  ice  grains,  evacuating 
and  flooding  the  mold  with  degassed  water  and 
then  freezing  the  resulting  ice-water  mixture. 
However,  details  of  the  procedure,  such  as  the 
source  of  the  seed  grains,  the  size  range  of  the  seed 
grains,  the  rate  and  direction  of  freezing  final 
specimen  porosity,  and  the  chemical  purity  of  the 
melt,  often  go  unreported.  These  factors  can  in¬ 
fluence  the  mechanical  behavior  of  the  material 
under  certain  conditions,  and  when  data  from  dif¬ 
ferent  sources  are  compared,  possible  specimen 
differences  must  always  be  considered.  In  addi¬ 
tion,  when  grain  size  is  varied,  questions  regarding 
the  influence  of  specimen  size  generally  arise  as 
well.  For  example,  Baker’s  (1978)  specimens  were 
19.7  mm  in  diameter  and  his  grain  diameter  ranged 
from  0.62  to  2.11  mm.  Duval  and  LeGac  (1980) 
used  specimens  of  80-mm  diameter  and  their  grain 
size  range  was  1.07  to  9.8  mm.  Thus,  although 
stress  and  temperature  were  the  same  for  both  sets 
of  data,  the  sample  sizes,  and  hence  the  number  of 
grains  across  a  sample  diameter  for  a  given  grain 
size,  vary  greatly. 

Unfortunately,  the  works  cited  above  are  not 
strictly  germane  to  the  present  study  because  of 
the  significant  difference  in  creep  stress  levels.  The 
higher  stress  level  generally  results  in  material  be¬ 
havior  in  the  dislocation  glide  with  cracking  re¬ 
gime  rather  than  a  diffusion  controlled  regime. 
Also,  the  grain  size  range  of  the  present  work  in¬ 
duces  a  significant  change  in  the  material’s  re¬ 
sponse  to  stress  by  causing  the  onset  of  internal 
fracturing.  Thus  the  grain  size  effect  found  here  is 
only  relevant  when  considering  the  onset  of  inter¬ 
nal  cracking.  A  study  of  grain  size  effects  near  the 
ductile-to-brittle  transition  offers  an  inherent  ad¬ 
vantage  since  the  effect  of  the  grain  size  variations 
is  evidenced  by  visible  cracking.  Thus,  a  grain  size 
effect,  and  the  associated  shift  in  deformational 
mechanism,  can  be  verified  visually.  This  is  a  pref¬ 
erable  situation  to  the  case  cited  above,  where 
stress-strain  rate  data  were  the  only  evidence  of  an 


apparent  deformational  mechanism  change,  and 
no  independent  means  of  verification  was  avail¬ 
able. 

In  Figure  16,  which  shows  minimum  creep  rate 
vs  the  average  grain  size,  there  appears  to  be  a 
subtle  trend  for  the  lowest  values  of  minimum 
creep  rate  to  occur  near  a  grain  size  of  3.0  mm. 
The  creep  rates  at  the  largest  grain  sizes  exhibit  a 
fairly  large  degree  of  scatter,  however,  making  it 
difficult  to  discern  a  trend  beyond  the  3-mm  grain 
size. 

It  is  interesting  to  note  that  the  more  clearly  de¬ 
fined  drop  in  f„,„  for  grain  sizes  between  1 .5  and  3 
mm  coincides  with  the  transition  from  the  thresh¬ 
old  of  cracking  to  a  significant  degree  of  cracking. 
Over  this  same  range  in  grain  size,  the  strain  at  f„,„ 
undergoes  a  significant  drop  from  10'2  to  approxi¬ 
mately  5.5xlO*3.  Thus,  the  material  appears  to 
lose  ductility  as  a  result  of  the  grain  size  increase. 
A  good  deal  more  testing  will  be  needed  to  clarify 
the  trends  in  these  creep  results. 

It  is  possible  that  the  behavior  seen  in  Figure  16 
is  merely  the  result  of  random  variations  in  the 
balance  between  competing  deformational  mecha¬ 
nisms.  This  is  an  inherent  problem  at  the  transi¬ 
tion  point  between  two  distinct  regimes  of  materi¬ 
al  behavior. 

Another  possibility  for  the  apparent  grain  size 
effect  for  the  larger  grain  sizes  relates  to  specimen 
size  effects.  When  grain  size  varies  as  in  the  pres¬ 
ent  work,  the  number  of  grains  in  a  specimen  of 
fixed  dimensions  varies  greatly,  as  does  the  num¬ 
ber  of  grains  across  the  diameter.  In  fact,  a  few  of 
the  larger-grained  specimens  tested  are  somewhat 
over  the  acceptable  limit  of  10  to  12  grains  across 
the  diameter.  (Note,  however,  that  when  using  the 
results  of  the  intercept  method  of  grain  size  esti¬ 
mation,  all  specimens  appear  to  be  within  the  limit 
of  10  grains  across  the  diameter.)  Jones  and  Chew 
(1983)  recommended  having  at  least  12  grains 
across  the  diameter  to  avoid  specimen  size  effects. 
Their  results  indicated  a  noticeable  increase  in  uni¬ 
axial  compressive  strength  when  the  number  of 
grains  across  the  specimen  dropped  to  eight.  In 
these  tests,  grain  size  was  held  constant  at  1 .0  mm 
and  the  specimen  size  was  changed  to  achieve  the 
range  in  the  number  of  grains  per  diameter.  There 
is  still  some  uncertainty  as  to  specimen  size  effects 
in  the  testing  of  ice  and  it  is  possible  that  the  pres¬ 
ent  testing  methods  do  not  completely  isolate 
grain  size  effects  from  possible  specimen  size  ef¬ 
fects.  A  considerable  amount  of  work  will  be 
needed,  however,  to  clarify  the  roles  of  grain  size 
and  specimen  size  in  the  mechanical  testing  of  ice. 
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Potential  difficulties  in  this  regard  are  recog¬ 
nized,  but  since  it  is  not  of  primary  concern  in  the 
present  work,  this  matter  will  not  be  dealt  with 
further. 

The  effect  of  grain  growth 

An  interesting  aspect  of  material  behavior  came 
to  light  regarding  the  effect  of  time/temperature 
history  on  the  creep  response  of  the  ice.  Specimen 
grain  size  is  generally  controlled  by  the  seed  grain 
size  and  a  given  sample  is  tested  soon  after  mold¬ 
ing  to  avoid  grain  growth  effects.  However,  a 
large  average  grain  size  can  be  achieved  by  suit¬ 
ably  aging  a  smaller  grain-sized  specimen.  The 
question  naturally  arises  as  to  whether  specimens 
of  equal  grain  size  display  similar  mechanical  be¬ 
havior  regardless  of  the  method  used  to  achieve 
the  grain  size.  Figure  15f  shows  data  which  ad¬ 
dress  this  point.  The  grain  sizes  of  specimens  41, 
49  and  56  in  Figure  1 5 f  were  achieved  by  allowing 
grain  growth  to  occur  for  some  time  after  mold¬ 
ing.  The  seed  grains  for  these  specimens  were  in 
the  0.59-0.83  mm  range.  The  grain  size  in  all  other 
tests  resulted  directly  from  the  seed  grain  size  and 
no  significant  grain  growth  occurred  before 
testing.  The  difference  in  the  creep  behavior  bet¬ 
ween  the  two  groups  of  specimens  is  striking. 
Specimens  41  and  49  do  not  exhibit  minimum 
strain  rates  as  such,  and  56  merely  develops  a  sub¬ 
tle  trend  near  10'!  strain  somewhat  indicative  of  a 
minimum  strain  rate.  The  common  trend  here  is 
the  absence  of  the  decreasing  strain  rate  usually 
found  in  primary  creep.  Additionally,  these 
specimens  always  show  an  extreme  degree  of  inter¬ 
nal  cracking. 

The  reasons  for  the  anomalous  behavior  of  the 
grain-growth  specimens  are  unclear.  It  is  unlikely 
that  a  preferred  orientation  developed  during 
grain  growth  under  these  conditions  when  the 
grains  were  originally  randomly  oriented.*  A  pos¬ 
sible  explanation  may  be  related  to  the  dislocation 
density  of  the  material  just  prior  to  testing.  The 
grain-growth  specimens,  after  being  aged  for  sev¬ 
eral  weeks  at  a  relatively  warm  temperature  (i.e. 
-2°C),  presumably  had  a  significantly  lower  dislo¬ 
cation  density  than  the  specimens  that  were  tested 
shortly  after  molding.  This  difference  in  disloca¬ 
tion  density  can  cause  a  corresponding  difference 
in  the  value  of  the  stress  needed  to  start  disloca¬ 
tion  motion  in  the  two  types  of  specimens.  Arm¬ 
strong  et  al.  (1962)  showed  that  the  frictional 
stress  term  an  increases  as  a  material  undergoes  the 
increase  in  dislocation  density  associated  with 
work  hardening.  If  this  is  the  case,  the  grain- 


Figure  27.  Creep  curve  for  a  fine-grained  specimen 
under  high  load  (a  =  4. 12  MPa). 


grov  th  specimens  will  experience  higher  effective 
shear  stresses  than  the  other  specimens.  This  in 
turn  leads  to  higher  internal  stress  concentrations 
and  hence  the  greater  degree  of  cracking.  Appar¬ 
ently,  this  greater  degree  of  damage  through 
cracking  is  associated  with  a  reduction  in  the 
strain  at  the  minimum  creep  rate  or,  in  some  cases, 
the  absence  of  a  discernible  minimum  creep  rate. 
It  is  expected  that  this  behavior  would  result  if  the 
applied  stress,  in  a  test  on  material  having  a  great¬ 
er  initial  dislocation  density,  was  sufficiently  high 
to  give  the  same  level  of  effective  shear  stress. 
Some  evidence  exists  in  support  of  this,  namely 
fine-grained  ice  without  any  grain  growth  was 
found  to  experience  only  a  brief  strain  rate  mini¬ 
mum  at  3.4  xlO'3  strain  under  a  stress  of  4.12 
MPa  (see  Fig.  27).  A  severe  amount  of  cracking 
accompanied  this  behavior.  Specimens  subjected 
to  slightly  lower  stresses,  however  (i.e.  3.7  MPa), 
exhibit  typical  creep  behavior  (see,  for  example, 
Mellor  and  Cole  1982)  and  strains  at  tm,„  are  near 
10  !.  Presumably,  the  strain  rate  minimum  would 
disappear  completely  under  some  further  increase 
in  stress.  Thus,  there  is  an  indication  that,  at  some 
level  of  effective  shear  stress,  ice  essentially  fails 
upon  loading  and  does  not  develop  the  strain  rate 
trends  typical  of  creep  behavior. 

An  additional  factor  complicating  an  assess¬ 
ment  of  the  observed  behavior  is  that  the  speci¬ 
mens  experiencing  grain  growth  are  likely  to  have 
a  rather  broad  range  in  grain  sizes  because  the 
larger  grains  grow  at  the  expense  of  the  smaller 
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a.  Thin  section  of  an  untested 
specimen  showing  an  extreme  ex¬ 
ample  of  the  grain  growth  pro¬ 
cess. 


b.  Thin  section  of  specimen  49 
after  testing.  Note  the  range  in 
grain  size  resulting  from  the  grain 
growth  process. 


Figure  28.  The  effect  of  grain  growth  on  grain  size. 


grains.  Figure  28a  shows  an  extreme  example  of 
this.  It  shows  an  untested  specimen  held  at  -2°C 
for  approximately  three  months  before  the  section 
was  taken.  It  exhibits  abnormal  grain  growth  as 
well  as  an  extremely  broad  range  in  apparent  grain 
diameters.  Figure  28b  shows  a  thin  section  of  spec¬ 
imen  49  after  testing  under  2.8  MPa  to  a  strain  of 
2.5  x  10~\  Note  that  some  grains  have  grown  con¬ 
siderably  while  clusters  of  fine  grains  (near  the 
center  of  the  photograph),  apparently  from  the 
original  structure,  still  persist.  It  is  difficult  to  de¬ 
termine  which  characteristics  of  such  a  structure 
control  the  deformational  processes. 

Due  to  the  uncertainties  involved  and  the  limit¬ 
ed  amount  of  data  available,  it  was  decided  not  to 
pursue  the  effect  of  grain  growth  on  mechanical 
behavior  in  the  present  work.  Once  the  above- 
mentioned  deviations  were  encountered,  speci¬ 
mens  were  tested  only  as  molded,  not  allowing  sig¬ 
nificant  grain  growth  to  occur. 

Normalized  crack  length 

It  is  useful  to  normalize  the  crack  length  data 
given  in  Appendix  A  to  the  grain  size  of  each  spec¬ 
imen.  This  allows  a  broad  comparison  of  the  re¬ 
sults  and  sheds  light  on  the  relationship  between 
the  crack  size  distribution  and  the  grain  size.  Fig¬ 
ure  29  shows  a  histogram  of  some  2246  observa¬ 
tions  made  from  the  thick  sections.  These  data 
have  been  normalized  to  grain  size.  The  mean  is 
0.5  and  the  standard  deviation  is  0.39.  Table  4 
gives  the  values  of  the  normalized  mean  crack 


length  for  all  specimens  (CL/d).  Although  these 
values  display  a  certain  amount  of  scatter,  they  are 
reasonably  well  grouped  about  the  mean.  Interest¬ 
ingly,  the  distribution  of  the  merged  normalized 
data  retains  essentially  the  same  shape  as  the  raw 
crack  size  data  of  the  individual  specimens  in  Ap¬ 
pendix  A.  Given  the  above  observations,  it  ap¬ 
pears  likely  that  a  generalized  distribution,  such  as 
that  in  Figure  29,  in  terms  of  normalized  crack 
length,  may  be  used  to  estimate  the  actual  crack 
size  distribution  for  any  given  grain  size. 

The  main  difficulty  in  this  connection  lies  in  as¬ 
sessing  the  maximum  crack  length.  There  is  con¬ 
siderable  scatter  in  the  largest  normalized  crack 
length  values  for  the  specimens  tested.  Figure  30 
shows  the  maximum  normalized  crack  length  as  a 
function  of  grain  size  for  all  tests.  They  range 
from  0.83  to  3.10  and  do  not  appear  to  correlate 
with  either  grain  size  or  axial  strain  level.  Perhaps 
fortuitously,  the  maximum  crack  size  of  3Ad 
agrees  with  the  observations  made  by  McMahon 
and  Cohen  (1965)  cited  earlier.  The  average  of  the 
maximum  normalized  crack  lengths  is  1.70  with  a 
standard  deviation  of  0.66.  Additionally,  the  raw 
data  indicate  that  8.8%  of  the  observed  cracks  are 
greater  in  length  than  the  average  grain  size. 

A  Beta-distribution  fit  to  all  the  data  indicates 
that  the  probability  of  encountering  a  crack  equal 
to  or  larger  than  1 .6 d  is  1  %  and  the  probability  of 
encountering  a  crack  equal  to  or  larger  than  2.0 d 
is  0.!°7o,  In  other  words,  it  is  relatively  rare  to  en¬ 
counter  a  crack  larger  than  the  average  maximum 
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Figure  29.  Normalized  fracture  length  distribution  for  alt 
tests. 
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Figure  30.  Maximum  normalized  crack 
length  vs  grain  diameter  for  all  tests. 


value.  In  fact,  all  observed  cracks  with  normalized 
values  greater  than  2.2  occurred  in  either  specimen 
71  or  74.  The  reason  for  the  unusually  high  values 
in  these  particular  specimens  is  not  apparent,  but 
the  likelihood  of  their  existence  must  nonetheless 
be  considered  in  applying  these  results. 

An  estimate  of  the  actual  crack  size  distribution 
of  material  with  an  arbitrary  grain  size  d  can  be 
obtained  from  the  normalized  crack  size  distribu¬ 
tion  by  substituting  the  term  CL/d  for  the  normal¬ 
ized  crack  size  CLN  and  then  multiplying  the  coef¬ 
ficient  by  l/d  to  maintain  unit  area  of  the  proba¬ 
bility  density  function. 

Location  of  cracks 

As  mentioned  above  (see  Table  3),  thin  sections 
of  several  highly  cracked  specimens  were  examined 
in  detail  in  order  to  assess  the  location  of  the  mi¬ 
crofractures.  The  microfractures  were  categorized 
as  either  grain  boundary  or  transcrystalline.  In 
total,  573  observations  were  made  and  Figure  31 


shows  the  results  in  the  form  of  histograms.  These 
data  have  been  normalized  to  grain  size. 

All  of  these  specimens  were  strained  to  10~J 
under  the  2.0-MPa  initial  creep  stress.  There  is  no 
apparent  systematic  variation  of  crack  location 
with  grain  size  under  these  conditions. 

The  mean  lengths  of  the  normalized  grain  boun¬ 
dary  and  transcrystalline  cracks  were  0.37  and 
0.35  respectively.  The  maximum  values  were  1.6 
and  1.2  respectively.  The  fact  that  these  measure¬ 
ments  came  from  thin  sections  probably  led  to  the 
lower  mean  and  extreme  values  relative  to  results 
obtained  from  thick  sections.  It  would  be  highly 
unlikely  for  the  maximum  crack  dimension  to  lie 
in  an  arbitrarily  selected  thin  section  of  1  mm, 
while  the  probability  of  observing  the  complete 
length  of  a  crack  in  a  thick  section  (of  10  mm)  is 
much  greater.  Thus,  the  absolute  magnitude  of  the 
results  in  Figure  31  should  be  treated  with  some 
caution. 
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a.  Normalized  grain  boundary  crack  size  distribution. 
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b.  Normalized  transcrystalline  crack  size 
distribution. 

Figure  31.  Normalized  crack  length  histograms. 


Acoustic  emission  activity 

In  the  present  work,  the  acoustic  emission  data 
serve  as  a  vehicle  to  link  the  fracturing  activity 
with  time  and  thus  specimen  strain.  The  initial 
task  concerns  the  development  of  a  correspond¬ 
ence  between  an  acoustic  event  and  the  nucleation 
of  a  discernible  microfracture.  The  correspond¬ 
ence  is  based  on  the  reasonable  assumption  that 
the  recorded  AE  event  amplitude  is  in  proportion 
to  the  microfracture  size,  as  noted  in  an  earlier 
section. 

The  crack  density  measurements  allow  the  esti¬ 
mation  of  the  total  number  of  visible  cracks  in  a 
given  specimen.  The  AE  monitoring  system  is  set 


to  a  sensitivity  great  enough  to  respond  to  acoustic 
activity  of  much  lower  amplitudes  than  that  gener¬ 
ated  by  the  observable  microfractures.  Thus,  in  a 
given  test,  there  are  generally  many  more  AE 
events  recorded  than  cracks  nucleated.  Conse¬ 
quently,  filtering  the  AE  data  was  necessary  in 
such  a  manner  as  to  retain  the  appropriate  number 
of  events  corresponding  to  the  estimated  total 
number  of  microfractures  in  the  specimen.  Events 
were  filtered  with  respect  to  amplitude  only.  A 
computer  program  performed  the  filtering  process 
in  two  modes:  1)  given  a  threshold  amplitude 
level,  it  determined  the  number  of  events  having 
amplitudes  greater  than  or  equal  to  the  threshold 
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using  a  simple  sorting  method,  and  2)  given  a  spe¬ 
cific  number  of  events,  it  determined  the  AE  am¬ 
plitude  threshold  that  was  passed  the  required 
number  of  times.  The  latter  mode  proved  most 
useful  in  the  present  context.  The  program  was  en¬ 
tered  with  the  estimated  number  of  microfractures 
and  an  output  file  was  in  turn  generated  that  con¬ 
tained  only  the  AE  events  that  passed  the  filtering 
process.  The  file  contained  the  time  of  occurrence 
and  the  amplitude  of  each  event. 

In  most  tests,  a  remotely  controlled  solenoid  im¬ 
parted  a  trigger  signal  to  the  specimen  at  the  mo¬ 
ment  the  load  application  began.  The  AE  system 
sensed  this  signal  and  its  time  of  occurrence  was 
taken  as  the  zero  or  reference  time  for  the  test. 
The  deformation  readings  also  were  “zeroed”  to 
this  time  to  ensure  a  common  starting  time  for 
both  the  AE  and  deformation  data.  Due  to  the 
method  of  recording  test  information,  two  sepa¬ 
rate  files  were  initially  developed  for  each  test. 
One  consisted  of  AE  data  as  a  function  of  time 
and  another  was  deformation  as  a  function  of 
time.  Both  files  were  interpolated  to  yield  readings 
at  the  same  time  increments  and  merged  to  form  a 
single  file. 

The  resulting  file  contained  sufficient  informa¬ 
tion  to  determine,  for  specific  time  increments,  the 
specimen  strain,  strain  rate,  stress,  accumulated 
fractures,  fracturing  activity  per  unit  time  and  per 
unit  strain.  The  acoustic  activity  was  normalized 
to  unit  volume  in  these  calculations. 

An  eventual  goal  of  the  AE  work  is  to  firmly  es¬ 
tablish  a  correspondence  between  acoustic  activity 
and  fracturing  activity.  This  relationship  will  al¬ 
low  a  prediction  of  the  size  distribution  of  the 
fractures  generating  the  observed  acoustic  activi¬ 
ty.  However,  at  this  writing,  this  relationship  is 
not  yet  sufficiently  established  to  warrant  a  de¬ 
tailed  analysis.  For  the  present  the  filtered  AE  re¬ 
sults  may  be  reasonably  taken  to  indicate  the  oc¬ 
currence  of  visible  cracks.  The  AE  observations 
for  each  test  are  based  on  filtered  data  employing 
the  filtering  thresholds  determined  for  each  speci¬ 
men. 

The  results  allow  the  determination  of  the  time 
and  strain  at  the  onset  of  fracturing  as  well  as  at 
the  peak  fracturing  rate,  as  given  in  an  earlier  sec¬ 
tion. 

With  the  exception  of  specimens  43,  49,  55  and 
70,  the  amplitude  thresholds  fall  within  the  range 
of  81 .6  to  88.7  dB,  with  a  mean  and  standard  devi¬ 
ation  of  86.2  dB  and  2.0  dB.  When  all  data  are 
considered,  the  mean  amplitude  threshold  is  83.7 
dB  with  a  standard  deviation  of  4.7  dB.  It  is  not 


clear  why  the  four  specimens  mentioned  have  sig¬ 
nificantly  lower  filtering  thresholds,  but  there  are 
several  sources  of  error  that  could  contribute  to 
inaccuracies  in  the  AE  measurements: 

1)  inconsistencies  in  the  characteristics  of  the 
specimen-transducer  interface 

2)  signal  attenuation 

3)  crack  orientation  effects 

4)  frequency  effects. 

A  variable  in  the  AE  considerations  for  all  tests 
is  the  consistency  or  repeatability  of  the  character¬ 
istics  of  the  specimen-transducer  interface.  Al¬ 
though  a  great  effort  was  made  to  be  as  consistent 
as  possible  in  the  placement  of  the  transducers, 
variations  in  the  quality  of  the  specimen/trans¬ 
ducer  interface  can  nonetheless  occur.  The  meas¬ 
ured  amplitude  of  identical  acoustic  pulses  de¬ 
creases  with  poor  specimen-transducer  contact. 
Ideally,  the  system  should  be  calibrated  with  an 
acoustic  pulse,  similar  to  that  generated  by  a 
nucleating  crack,  just  prior  to  each  test.  Thus,  de¬ 
ficiencies  could  be  detected  beforehand  and  the 
transducer  remounted  to  provide  a  satisfactory  re¬ 
sult.  Unfortunately,  such  a  method  was  not  avail¬ 
able  during  the  course  of  this  study,  so  some  un¬ 
certainty  is  inherent  in  the  AE  results. 

In  fact,  such  uncertainty  could  be  the  major 
cause  of  variations  in  the  filtering  thresholds  (see 
Table  5).  This  circumstance  causes  difficulty  in  the 
longer-range  objective  of  establishing  an  overall 
amplitude  threshold  for  visible  fractures.  How¬ 
ever,  it  does  not  adversely  affect  the  veracity  of 
the  results  when  specimens  are  considered  individ¬ 
ually.  Thus,  even  though  the  thresholds  may  vary, 
the  filtered  AE  events  for  a  particular  specimen 
correspond  to  the  observed  number  of  cracks  for 
that  specimen.  Consequently,  such  quantities  as 
the  cracking  rate  and  the  strain  level  for  the  onset 
or  peak  rate  of  cracking  are  not  directly  influenced 
by  small  variations  in  the  specimen-transducer  in¬ 
terface  quality. 

Attenuation  of  the  acoustic  signal  is  significant 
for  high  frequencies  and  large  distances  in  ice. 
However,  in  this  study,  the  maximum  travel  dis¬ 
tance  from  an  event  source  to  one  of  the  trans¬ 
ducers  is  approximately  50  mm.  Bogorodskii  and 
Gusev  (1973)  indicate  that  attenuation  is  on  the 
order  of  5  dB  m'1  in  ice.  Thus,  signal  variations  of 
much  less  than  1  dB  are  expected  in  the  present 
work,  and  are  thus  not  a  significant  factor  in  the 
results. 

The  orientation  of  the  crack  to  the  transducer 
face  is  likely  to  have  an  effect  on  the  measured 
amplitude.  The  effect  is,  however,  difficult  to 
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Figure  32.  Mean  AF  amplitude  vs  mean  crack  length. 


quantify.  It  is  assumed  that  the  orientation  of  the 
cracks  to  the  transducer  face  is  sufficiently  ran¬ 
dom  to  preclude  systematic  effects  of  this  varia¬ 
ble. 

This  treatment  does  not  address  the  effects  of 
possible  variations  in  the  frequency  content  of  the 
acoustic  pulses.  Since  the  mechanics  of  crack  for¬ 
mation  are  expected  to  be  the  same  for  all  crack 
sizes,  little  variation  is  expected.  However,  the 
possibility  that  the  frequency  spectra  of  the  acous¬ 
tic  pulses  may  change  with  crack  size  should  be 
kept  in  mind. 

As  pointed  out  in  an  earlier  section,  a  micro- 
fracture  is  expected  to  generate  an  acoustic  pulse 
in  proportion  to  its  size,  with  certain  geometric 
considerations.  Thus,  the  magnitude  of  the  acous¬ 
tic  pulses  passing  the  filter  should  correspond  to 
the  sizes  of  the  observed  microfractures.  Figure  32 
shows  the  average  AE  amplitude  vs  the  average 
crack  size  for  all  the  specimens  for  which  suitable 
AE  data  were  obtained.  Included  in  this  plot  is  a 
point  from  specimen  69,  which  evidenced  no  visi¬ 
ble  cracking.  This  amplitude  of  79  dB  was  very 
rarely  exceeded  in  the  test  and  is  assumed  to  repre¬ 
sent  the  greatest  AE  amplitude  produced  by  a  sub- 
visible  crack. 

These  data  indicate  the  tendency  for  the  mean 
AE  amplitude  to  increase  with  mean  crack  size. 
The  scatter,  however,  is  considerable  and,  as  men¬ 
tioned  above,  this  scatter  is  probably  due  to  varia¬ 
tions  in  the  quality  of  the  specimen/transducer  in¬ 
terface.  The  average  crack  size  -  average  ampli¬ 
tude  relationship  could  be  strengthened  by  an  in¬ 
crease  in  the  relatively  narrow  range  of  0.90  to  2.5 
mm  in  the  average  crack  size.  Significantly  larger 
specimens  would  be  required  to  increase  this  range. 

The  onset  of  cracking,  as  indicated  by  the  AE 
results,  occurs  at  an  axial  strain  of  4.7  x  10  4  with  a 
standard  deviation  of  4.3  x  10'4.  Interestingly,  this 


strain  level  is  in  good  agreement  with  that  given  by 
Gold  (1967)  for  the  strain  at  the  formation  of 
grain-sized  cracks  in  columnar-grained  ice.  The 
strain  and  time  at  the  onset  of  cracking  activity  do 
not  appear  to  vary  systematically  with  grain  size  in 
these  tests.  Strain  values  for  the  cracking  onset 
range  from  effectively  zero  to  a  maximum  of 
1.15x10'. 

As  Figure  23  shows,  grain  size  exerts  a  strong  in¬ 
fluence  on  the  time  to  the  maximum  cracking  rate 
(as  indicated  by  the  maximum  AE  rate).  Since  the 
strain  at  \dAE/dt\„,„  is  not  strongly  affected  by 
the  grain  size  in  these  tests,  it  is  evident  that  the 
average  strain  rate  prior  to  the  peak  cracking  rate 
must  increase  with  grain  size.  Indeed,  an  inspec¬ 
tion  of  the  creep  curves  (Fig.  15c)  shows,  for  a 
constant  initial  stress  of  2.0  MPa,  that  the  average 
strain  rate  prior  to  [dAE/di]„ax  increases  by  nearly 
a  factor  of  four  as  grain  size  increases  from  1 .8  to 
5.5  mm.  There  is  a  very  moderate  tendency  for  the 
strain  at  \dAE/dt]mxx  to  decrease  with  increasing 
grain  size,  further  contributing  to  observed  de¬ 
crease  in  time  to  the  maximum  cracking  rate.  The 
average  strain  at  [dAE/de]ma.  is  1.95x10'’  with  a 
standard  deviation  of  9.6  xlO  4.  The  average 
strain  at  [ dAE/dt\ma ,  is  1 .75  x  10'’  with  a  standard 
deviation  of  8.8  x  10\ 

In  general,  these  results  reinforce  earlier  find¬ 
ings  regarding  cracking  in  polycrystalline  ice.  A 
small  amount  of  strain  is  required  to  initiate 
cracking,  and  once  started,  the  rate  of  cracking 
reaches  a  peak  during  the  primary  creep  stage, 
well  before  the  minimum  creep  rate  is  reached. 

The  average  strain  at  which  [dAE/dt]m.x  occurs 
coincides  with  the  average  strain  at  the  in¬ 
flection  point  found  in  the  primary  creep  portion 
of  a  plot  of  e  vs  f  (see  Fig.  15).  (Both  these  values 
of  average  strain  are  very  near  2x10"’.)  In  the 
present  deformational  mechanism  regime  of  plas- 
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tic  flow  with  cracking,  AE  results  generally  show 
no  distinct  characteristic  at  either  the  occurrence 
of  in  a  creep  test  or  in  a  strength  test  (see 
Cole  and  St.  Lawrence  1984).  However,  it  is  ex¬ 
pected  that  the  maximum  cracking  rate  should 
coincide  with  a  fundamentally  significant  aspect 
of  material  behavior.  Mellor  and  Cole  (1982)  were 
the  first  to  distinguish  the  point  of  in  ice  creep 
data,  and  noted  that  it  is  the  point  of  maxi¬ 
mum  deceleration  in  the  primary  creep  curve. 
Thus,  the  AE  results  indicate  that  the  material  is 
apparently  undergoing  its  greatest  rate  of  strain 
hardening  (i.e.  creep  rate  is  decreasing  most  rapid¬ 
ly)- 

Although  this  application  of  AE  technology  re¬ 
quires  further  development,  the  potential  for  its 
use  is  clear.  Proper  interpretation  of  AE  data  will 
preclude  the  time-consuming  post-test  analysis 
currently  required  for  examining  internal  cracking 
activity.  Much  insight  can  be  gained  from  detailed 
information  on  the  characteristics  of  internal  frac¬ 
tures  in  ice  in  terms  of  understanding  deformation 
mechanisms  and  for  verification  of  micromechan¬ 
ical  models  of  material  behavior. 


SUMMARY  AND  CONCLUSIONS 

This  work  presents  the  results  of  constant  load 
creep  tests  performed  at  -5°C  on  polycrystalline 
ice.  Some  26  tests  were  performed  on  specimens 
having  equiaxed  grains  ranging  from  1.5  to  6.0 
mm.  Most  specimens  experienced  an  initial  stress 
of  2.0  MPa.  Tests  were  terminated  after  axial 
strains  ranging  from  3.7x10  4  to  5x10*.  Some 
specimens  were  tested  at  higher  stress  levels  of  2.4, 
2.6  and  2.8  MPa. 

The  results  demonstrate  the  influence  of  grain 
size  on  the  internal  fracturing  in  ice.  The  stated  in¬ 
crease  in  grain  size  brought  about  the  onset  of  in¬ 
ternal  cracking  as  predicted  by  the  nucleation 
equation  of  Smith  and  Barnby  (1967).  The  materi¬ 
al  showed  a  loss  in  ductility  as  evidenced  by  a  sig¬ 
nificant  decrease  in  the  strain  at  and  a  dram¬ 
atic  increase  in  the  number  of  internal  cracks. 

In  all  tests  where  internal  cracking  occurred,  an 
extensive  post-test  analysis  on  the  size  and  number 
of  cracks  in  the  ice  was  made.  This  analysis  al¬ 
lowed  the  determination  of  the  crack  size  distribu¬ 
tion  as  well  as  the  crack  density  in  each  specimen. 
A  linear  relationship  between  grain  size  and  the 
average  crack  size  distribution  emerged  from  these 
data. 

Peripheral  aspects  of  the  work  addressed  crack 
healing,  the  theoretical  aspects  of  the  grain  size/ 


crack  size  relationship,  the  grain  size  effect  on  the 
strain  at  <«,„  the  relationship  between  grain  size 
and  slip  plane  length  and  the  effect  of  grain  size  on 
creep  behavior. 

Acoustic  emissions  monitoring  techniques  were 
employed  and  the  results  were  promising.  Consid¬ 
erable  progress  was  made  toward  determining  the 
onset,  total  number  and  rate  of  formation  of  mi¬ 
crocracks  from  AE  data. 

In  view  of  the  test  conditions  stated  above,  this 
work  leads  to  the  following  conclusions: 

1.  Both  the  crack  nucleation  condition  and  the 
crack  size/grain  size  relationship  are  well  modeled 
using  the  concept  of  the  dislocation  pileup. 

2.  Cracks  begin  to  nucleate  in  ice  at  grain  sizes 
between  1.5  and  1.8  mm  under  a  2.0-MPa  creep 
stress  at  -5  °C. 

3.  The  extent  of  internal  microfracturing  in¬ 
creases  sharply  as  grain  size  increases  from  1.8  to 
6.0  mm. 

4.  The  average  microcrack  dimension  scales 
linearly  with  the  mean  grain  size. 

5.  As  grain  size  increases  from  1.5  to  6.0  mm, 
the  axial  strain  at  the  minimum  creep  rate  falls 
from  10'!  to  near  4x10"*. 

6.  The  minimum  creep  rate,  while  exhibiting 
some  subtle  trends,  is  not  clearly  affected  by  grain 
size. 

7.  The  peak  fracturing  rate  occurs  at  relatively 
low  strains,  well  before  tm,„  is  reached,  and  the 
time  to  the  peak  fracturing  rate  decreases  with  in¬ 
creasing  grain  size. 

8.  Grain  size  apparently  does  not  systematical¬ 
ly  affect  the  strain  at  the  onset  of  cracking  or  the 
strain  at  the  peak  fracturing  rate. 

9.  Based  on  a  limited  number  of  observations, 
grains  that  develop  transgranular  cracks  do  not 
appear  to  have  a  preferred  orientation  to  the  stress 
axis. 

10.  Crack  healing  processes  significantly  affect 
crack  size.  Water-vapor-filled  cracks  heal  by  a  vis¬ 
cous  flow  mechanism  at  a  significantly  higher  rate 
than  air-filled  cracks. 

11.  The  crack  lengths  for  all  specimens,  when 
normalized  to  grain  size,  tend  to  follow  a  common 
distribution. 


SUGGESTIONS  FOR  FUTURE  WORK 

1.  Experimental  work  must  be  carried  out  to  ex¬ 
plore  specimen  size  effects  on  the  mechanical  be¬ 
havior  of  polycrystalline  ice.  This  should  be  done 
in  conjunction  with  a  grain  size  study  in  order  to 
isolate  grain  size  and  specimen  size  effects. 
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2.  Both  strength  and  creep  tests  on  specimens  of 
varying  grain  size  should  be  performed  both  above 
and  below  the  threshold  of  cracking  as  indicated 
in  this  work.  This  would  help  clarify  the  trends 
suggested  in  the  present  data.  Specifically,  such  re¬ 
sults  would  show  if  grain  size  exerts  an  influence 
on  the  stress  vs  strain  rate  relationship  to  accom¬ 
pany  its  influence  on  the  internal  fracturing  activi¬ 
ty- 

3.  Continued  efforts  toward  establishing  a  crack 
size  vs  AE  amplitude  relationship  appear  warrant¬ 
ed.  A  useful  relationship  of  this  type  would  greatly 
reduce  the  amount  of  work  required  to  analyze 
fracturing  activity.  The  optical  post-test  analysis 
methods  used  in  this  work  are  too  time-consuming 
and  tedious  to  become  part  of  a  routine  testing 
procedure.  However,  the  AE  technique,  once  suf¬ 
ficiently  advanced,  will  lend  itself  to  routine  use. 

4.  A  detailed  study  of  crack  orientation  to  the 
axis  of  applied  stress,  especially  under  uniaxial 
tension,  would  be  useful.  Given  a  sufficient  num¬ 
ber  of  observations,  it  is  possible  to  establish  a 
distribution  of  the  angles  between  the  crack  plane 
and  the  stress  axis.  Knowledge  of  this  distribution 
as  well  as  its  development  during  the  course  of 
straining  will  shed  light  on  the  failure  process  of 
ice  in  tension. 

5.  Cracking  in  specimens  having  grains  signifi¬ 
cantly  larger  than  those  tested  in  this  work  should 
be  examined.  It  is  possible  that  subgrain  structure 
might  become  a  factor  at  large  grain  sizes  and  thus 
limit  the  applicability  of  the  crack  size  vs  grain  size 
relationships  found  in  the  present  work. 

6.  The  effects  of  grain  growth  on  the  mechani¬ 
cal  properties  of  ice  should  be  examined  system¬ 
atically.  It  is  possible  the  time-temperature  history 
should  be  considered  as  a  test  variable  along  with 
structural  parameters. 
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APPENDIX  B:  CRYSTAL  ORIENTATIONS 


This  appendix  contains  pole  diagrams  for  two  specimens  containing  cracks.  For  each 
specimen,  a  number  of  measurements  were  taken  to  indicate  the  random  nature  of  the  grain 
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